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ABSTRACT 


An expanded package of interactive FORTRAN computer pro- 
grams has been developed for the analysis and forecasting 
of time series data. The package, called the Time Series 
Editor, is designed to employ the iterative Box-Jenkins 
methodology of time series analysis. The Time Series Editor 
was developed for time-shared use on the Control Program/ 
Cambridge Monitor System (CP/CMS) at the U.S. Naval Post- 
graduate School, but can be modified for use on other time- 
sharing systems with a FORTRAN capability. The Time Series 
Editor assists in data preparation and entry, analysis, 
modeling, forecasting and diagnostic testing. Utilization 
of the package, following the included User's Guide, requires 
only a limited knowledge of the computer system, with all 


required user responses interactively prompted by the Editor. 
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Ley = INTRODUCTION 


Operations researchers, military planners and programmers, 
Statisticians, Beonemaetee marketing personnel, managers and 
@thers are often faced with the need to analyze data in the 
form of a time series, which can be thought of as a sequence 
Gt observations Of either a deterministic or stochastic 
BEOCESS - In most cases, the objective of the analysis is 
to determine patterns or other recognizable behavior appar- 
ent in the data, and to then formulate a suitable mathemati- 
cal model for the time series from which forecasts of future 
behavior can be obtained. The value to a decision maker 
of being able to predict the future with some reasonable 
and statistically quantifiable degree of accuracy can not 
be overstated. For example, in areas such as budget expen- 
ditures, recruiting performance, commodity prices, population 
levels, resource consumption, manpower levels and consumer 
demand for products, decision makers charged with planning 
for the future should base their decisions at. least in 
part on the best available predictions about the future 
behavior of the time series in question. 

Until the late 1960's, the techniques employed in time 
series analysis were primarily those of spectral analysis, 
with applications of harmonic analysis and mathematical 
EGanstorm theory. For an introduction to these topics, the 
reader is referred to Anderson [Ref. 2], Lewis [Ref. 7] 


and Jenkins [Ref. 10]. Due to the high degree of mathematical 
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sophistication required by the spectral analysis approach, 
the capability to perform time series analysis resided nearly 
exclusively with mathematigians and electrical engineers. 
As a consequence, the majority of decision makers came to 
use more understandable (and far less powerful) methods 
such as simple moving averages or exponential smoothing. 
However, in the late 1960's and early 1970's, the statistical 
analysis of time series by the methods developed by Box 
and Jenkins [Ref. 4] has gained widespread acceptance. 
Although these methods are mathematically nearly equivalent 
to the spectral approach (there are transformations that 
interconnect the two methodologies), the Box-Jenkins approach 
is described in vocabulary more familiar to operations 
researchers, statisticians, economists and managers, and 
is therefore being used more and more by them in building 
models from measurements of the past. 

Many algorithms and computer programs for performing 
the analyses required by the Box-Jenkins approach have been 
developed and are available from several sources. Perhaps 
the best source is the collection of FORTRAN computer sub- 
routines which resides in the International Mathematical 
and Statistical Library (IMSL) [Ref. 9]. The primary prob- 
lem with using the available computer resources lies not 
in any deficiency of the programs themselves, but with the 
very nature of the Box-Jenkins methodology itself. The 


Box-Jenkins method is an iterative approach to time series 
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analysis, which is briefly described in Figure l. [See 
Wheelwright and Makridakis, Ref. 16.] 

As indicated by Figure 1, the Box-Jenkins method is a 
multi-stage, iterative process. It begins with the postu- 
lation of a general class of time series models which has 
been found, experimentally, to be extremely rich. There- 
after, the procedure continues as a trial-and-error process, 
with decisions points where the analyst is required to 
select the next direction based on the best information 
available to him. Since each stage of the process outlined 
in Figure 1 may consist of several sub-steps, the modeling 
process itself can become quite time-consuming, even with 
the ready availability of the IMSL software resources. For 
example, an analyst employing the IMSL subroutines in a 
batch processing computer system to perform Box-Jenkins 
modeling of a time series might perform the following 
sequence of tasks: 

1. Prepare the time series data in the proper format. 

2. Plot and visually examine the time series, checking 
for nonstationarity, trends, seasonality, patterns, 
etc. 

3. Write a program to call the IMSL subroutine that 
will calcuate the mean, variance, autocorrelations 
and partial autocorrelations of the series. 

4. Plot the autocorrelations and partial autocorrela- 
tions; this provides much of the information re- 
quired for identifying the correct model for the 


series. 
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5. Write a program to call the IMSL subroutine which 
transforms the time series to adjust for seasonal 
patterns, nonstationary behavior or other behavior 
which deviates from that assumed for the class of 
models postulated. 

6. Repeat steps (2) through (4) with the transformed 
data. 

7. Review the statistical properties of the auto- 
correlations and partial autocorrelations for 
tentative model identification. 

8. Write a program to call the IMSL program that 
estimates the model parameters and computes the 
residuals. 

9. Write a program to perform goodness-of-fit tests 
for the model. 

10. Analyze the model residuals, using steps (1) through 
(9), as with the original time series. 

ll. Refine the model using information gained by 
examination for any structure in the residuals. 

12. Repeat the preceding process until no structure 
remains in the residuals. 

13. When an adequate model has been obtained, write a 
program to call the IMSL subroutine which forecasts 
and determines confidence intervals for future 
values of the time series. 

Between each succesSive pair of steps, the analyst must 


Manually intervene in the analysis process, making a decision 
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subjectively based on the information then available. 

Thus, much analyst interaction is necessary in order to 
determine a Suitable mathematical model and forecast equa- 
tion. Even with rapid computer job turnaround time, the 
process described above could easily consume a full working 
day or more. 

This paper describes an effort to eile aa aee some of the 
Organizational problems inherent in Ener aisaaduine of time 
series using the Box-Jenkins iterative approach. This 
effort consists of an interactive computer program package 
which provides access in a structured way to the most useful 
IMSL and other subroutines necessary for Box-Jenkins time 
series analysis. This package, called the Time Series Editor, 
is written for use with the Naval Postgraduate School's 
Control Progam/Cambridge Monitor System (CP/CMS). Since all 
working programs with the exception of the executive routine 
itself are written in FORTRAN, the Time Series Editor should 
be adaptable for use on other FORTRAN-capable time-sharing 
systems. The Time Series Editor assists the analyst in 
data preparation and entry, model construction and diag- 
nostic testing, and time series forecasting. In fact, 
with the User's Guide provided as an Appendix to this report, 
a complete Box-Jenkins time series analysis can be per- 
formed in a short time by even a novice computer user. 

For this reason, the Time Series Editor could be a valuable 
instructional tool for laboratory use in a time series 


course. 
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A brief overivew of Box-Jenkins methodology is provided 
in Chapter II, to serve as a point of reference for the 
material which follows. Chapter III contains a description 
of the algorithm employed to estimate non-linear least- 
squares parameters for a generalized (seasonal, nonseasonal, 
stationary, or nonstationary) Box-Jenkins model. Chapter IV 
contains descriptions of each of the major programs in 
the Time Series Editor. Chapters V and VI contain examples 
of actual use of the Time Series Editor, using a non- 
seasonal and a seasonal series, respectively. Chapter VII 
includes a summary of the report, and recommendations for 
future additions to the Time Series Editor. The Appendices 
include a User's Guide to the Time Series Editor, with 
instructions, sample user sessions and sample outputs, as 
well as listings of the Editor itself and all primary 


subprogram packages. 
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II. BOX-JENKINS METHODOLOGY: A PRIMER 


This chapter provides a brief overview of the methodology 
developed for the analysis of time series data by G.E.P. 
Box and G.M. Jenkins, in order to facilitate understanding 
of the models and programs employed in the Time Series 
Editor. For more detail on the material covered in this 
chapter, the interested reader is referred to the texts by 
Anderson [Ref. 1], Box and Jenkins [Ref. 4], Pindyck and 


Rubinfeld [Ref. 15], Nelson [Ref. 14], and Mabert [Ref. 12]. 


A. LISTING OF NOTATION 
The following is a listing of the notation employed 
in this chapter and in the remainder of this report. Each 


symbol is provided with a brief title or explanation of its 


use. 
{z,} a time series 
N the length of the time series (number of terms) 
Z, 62 estimates of the mean and variance of the time 
7 series {z, } 
e 
: r 
' —— + 
{zi} transform of {zt}, where z, (2, wc) 
’ = + 
or 2) In(z, G) 
iw} differenced series from {z,}, where 
wo = voyez. 
ce Sade 
P. aucocorrelation of lag k 
rx estimate of autocorrelation at lag k 
vt backward difference operator of order d 
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seasonal backward difference operator of order 
D with season s 


backward shift operator of order s 
length of the time series' season 


non-seasonal autoregressive parameters 
(estimates, 6) 


seasonal autoregressive parameters (estimates, ¢) 


non-seasonal moving average parameters 
(estimates, 96) 


seasonal moving average parameters (estimates, 0) 
series of white noise terms of the process 
variance of the white noise process 

om of squared residuals of the model 

number of non-seasonal autoregressive parameters 
number of seasonal autoregressive parameters 
number of non-seasonal moving average parameters 
number of seasonal moving average parameters 
transfer function of a linear filter 


partial autocorrelation of log k 


B. TIME SERIES, STOCHASTIC MODELS AND STATIONARITY 


A distinguishing feature of modern time series analysis 


is that the sequence of observations of a given variable is 


considered to be a realization of jointly distributed 
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random variables, such that the time series itself can be 
viewed as a stochastic process. A discrete stochastic 


time series, then, is a set of observations (2) 1Z5r-+-72,) 


generated sequentially in time by a set of jointly distributed 
random variables, where the actual data (Zy1Zor---7Z,) 
represents a particular realization of some joint probability 


distribution function E(Z)1Zor-2+ 72 ed With this function 


e 


determined in some way, a forecast Zz can be thought of as 


Ek 


having been generated by a conditional probability distri- 
bution £(Z, |Z, 02g0e++ 024) Now, the stochastic process 
which generates the time series (ZyrZor---+72,) is said to 

be stationary if its properties are unaffected by a change 

in time origin; that is, the process has reached a particular 
meate Of Statistical equilibrium. This implies that the ° 
joint probability distribution associated with the m 


observations (z ), Made at any set of times 


ida eR ay 
S| © en 
(tyrtgr---rt)e is identical to that distribution associated 
with the m observations (z ) made at 


12 yee 
t,+k t,tk 
SRS Pee neip te +k). 

m 


Tee ete 
mn 


times (t,+k,t 


a 2 
For stationary discrete time series, the marginal 
probability distribution p(z,) is the same for all times 


t. Hence, the stochastic process embodied in the time 


series has a constant mean 
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which defines the level about which it fluctuates in time, 


and a constant variance 


co 


ete a=)  - )-piz)az , 


which measures its variability about the level of the mean. 
Since the probability distribution ae the same for all 
times t, its form can be discovered by plotting a histogram 
or relative frequency plot of the set of observations 

(Zp rZor-++rZ,)- Additionally, the mean and variance of 

the stochastic process can be estimated by the sample 


averages taken over time: 


N 
— oe 
Z = N Daze 
=1 
and 
N 
AS pete —,2 
Soe y (22i5 "Zar: 
t=1 


C. THE BASICS OF DIFFERENCE OPERATORS AND LINEAR FILTER 
MODELS 


mais Section provides a brief discussion of the differ- 
ence operators currently employed in writing Box-Jenkins 


models. Those most commonly used are described below in 


Table I. 
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Operator Name Definition 


m 


B Backward Shift Operator BZa = 2 
tr t=1 
Of order m e 
B Z = Zeke 
Fo Forward Shift Operator F = Bee 
of order m 
P2e Rectan 
m 
FE Z = cara 
Ve Backward Difference z = (1 - BD 
Operator of order D 
and season s ve outa B2)? 
V2, Cer ier oa) 


Table I. Basic difference operators. 


It is noted here that ae could be written as a 


binomial expansion; that is, 


D 

D Sees a 
(eo (3) (2) (-B) eas 
j=0 


The stochastic models employed in the Box-Jenkins method 


are based on the idea that a time series in which successive 
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values are highly dependent can be regarded as having been 
generated from a series of independent "shocks" ay: These 
shock terms, sometimes referred to as white noise terms, 
are random drawings from some fixed distribution, generally 
assumed Normal with mean zero and variance ee The white 


noise process a, is supposed transformed to the process z 


& t 


by what is known as a linear filter, shown in Figure 2 below 


[Ref. 4, p. 8]. The linear filtering operation takes a 
v (B) 
White a es Process 
Noise Pilte: Zo 
fe 
Figure 2. Representation of a time series as 


the output from a linear, filter. 


weighted sum of previous observations, such that Zh can be 


written 


ll 
i= 

+ 
< 
tw 
ny) 


In general, wu is the parameter that determines the "level" 


of the process, and 
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y(B) = 1+y,B+ vB panes 


is the linear operator that transforms a, into z called 


t ce 
the transfer function of the filter. The sequence 
Vy Yas tHe formed by the weights may be (theoretically) 
finite or infinite; if this sequence is finite, or infinite 
guidmeonverges, ‘thesfiltersis said to be stable and the 
process Zz, Stationary. The parameter u is then the value 
about which the process varies. Otherwise, Zy is called 


nonstationary, and un has no specific meaning except as a 


reference point for the level of the process. 


D. HOMOGENEOUS NON-STATIONARY AND SEASONAL PROCESSES 

In many cases of interest to the time series analyst, 
the process being examined is not stationary. Instead, 
the probabilistic structure of the process generating the 
time series changes with time. For example, there may be 
a trend or seasonal pattern in the time series. However, if 
the series can be found to exhibit behavior which is somewhat 
homogeneous, then the series can frequently be transformed 
into a stationary series that can be described by Box-Jenkins 
models. 


A series is said to be homogeneous nonstationary of 


ender d if Tike ee ae is a stationary series for some integer 
d. Here V = (1 - B) denotes the difference operator such 
that Via, = 2 — Z2,_,- The fact that a given series is 


nonstationary can be determined quickly by examination of 
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the autocorrelation function, to be discussed later. The 
series to be modeled should be differenced until the re- 
sulting series appears stationary or until the procedure 
appears to be making no improvement. 

A time series is called seasonal when it exhibits 
cyclical behavior over time, showing a regular period. 
Examples of seasonal processes include monthly rainfall, 
monthly crop yields, sale volume of bathing suits, livestock 
production rates and energy consumption. Seasonal patterns 
are often easy to spot simply by observing a plot of the 
series, or through knowledge of the process that generated 
the series. However, in many cases where the variability of 
the series is large, seasonal patterns will be difficult to 
distinguish from other fluctuations. Recognition of season- 
ality in a series is important, since it provides information 
useful in modeling and forecasting; as will be shown, the 
autocorrelation function can assist in recognizing seasonality. 
The Box-Jenkins modeling approach for seasonal nonstationary 
time series of season length s is to first transform the 
series using an appropriate seasonal differencing operator 


to "Sweep out" the seasonal effect; that is, take 


If there is a trend or other types of non-stationarity pre- 


sent, as well as the seasonal effect, the series can be 
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differenced again with the ordinary difference operator 
until stationarity is achieved. This transform would be 


written 


The transformed series Wy is then modeled as a stationary 
Series. For seasonal series where complete stationarity 
can not be achieved, a class of seasonal models is avail- 


able; these are discussed later in this Chapter. 


E. GENERAL CLASSES OF BOX-JENKINS MODELS 

This section will provide brief descriptions of several 
of the most commonly used classes of Box-Jenkins models, 
including autoregressive (AR) models, moving average (MA) 
models, mixed and integrated autoregressive-moving average 
(ARMA/ARIMA) models, and a generalized seasonal ARIMA model. 

1. Autoregressive Models 

In autoregressive (AR) models, the current value of 

the series under study is expressed as a linear combination 
of previous series values that explain the current obser- 
vation, plus an unexplained random (white noise) term a,- 
For a stationary series, or one that has been transformed 
to stationarity, the deviation from the mean value for each 
period, 2, = Z, - uy, can be modeled as dependent on weighted 


values of the previous deviations from the mean. The AR 


model can be written 
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PREC eal ae eto? 
In the expressions which follow the tilde will be dropped, 
but it is to be understood thatthe mean has already been 
Subtracted from each observation. This equation is called 
an autoregressive process of order p, where ve represents 


h 


a scalar weighting coefficient for the 3° previous period. 


The model can also be written in more compact form as 


>(B) 2, Sees 


where 
7 2 Bi, s a! p 
o(B) = (1 18 >5B Beare 8 ae 


a polynomial in the operator B. Then, $(B) represents a 
general non-seasonal AR operator. The reason these models 
are called autoregressive can be seen from recalling 
regression analysis; here the model relates the dependent 
variable zy to a set of explanatory variables (2,372, 97 


Melep 2 ), which are previous values of the time series. 


t-p 
Therefore, the model is autoregressive. For example, when 
the value of the time series at time t depends only on the 
value at time t-l, the model becomes an autoregressive model 


Seporder 1, denoted AR(1). Similarly, if the value in 


period t depends on the values in periods t-l, t-2 and t-3, 
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2. Moving Average Models 


In moving average (MA) models, the assumption is made 


that the current value of the time series at period t can 
be expressed as a linear combination of the previous fore- 
cast errors (or, residuals), a, = Zy - Ze - The general 


equation for this model can be written 


which is called a moving average model of order q, where 
a represents the scalar weighting coefficient for the a 
previous period. Like the AR model, the MA model can be 


written in more compact operator notation as follows: 


where 
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a polynomial in the operator B. The MA model implies that 
the analyst can gain valuable information for Luture 
predictions by consideration of the weighted sum of a 
number of previous forecast errors or residuals. The 
following are examples of MA models of order 1 and order 


3, MA(1) and MA(3): 


MA(1) : Zz, o= (eS 8,B)a, = ae 05a.) ; 
MA(3) : SiCleceE seas scan- a 
bea I? 3 t 
a oeee Vere me Goce ee ua eR—3 


3. Mixed and Integrated Autoregressive-Moving Average 
Models 


A natural extension of the autoregressive and moving 


average models is to construct a combination of the two. 


Such mixed processes are called autoregressive-moving average 


(ARMA) processes, of orders p and q, often written as 
ARMA(p,q). Box and Jenkins [Ref. 4] notes that for many 
series encountered in practice, the inclusion of both AR 
and MA terms in a model results in fewer total parameters 
than would be required for either a pure AR or pure MA 
process. The ARMA(p,q) model may be written in the usual 


form as 
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in the more compact operator notation, the same model is 


written 
o(B)Zz, = 6(B)a 
As an example, an ARMA(2,2) model would be written as: 
Pi hr sete) O20 hg wel ees Bye 
it 2 t 
or expanded, as: 


Reet eee 3-E=2 oe ka eH? 
Recall that with differencing or other transformation 
techniques, a stationary series can be obtained from a 
homogeneous (Seasonal or non-seasonal) non-stationary series. 
When the differencing technique is integrated directly into 
the ARMA model, the result is an integrated autoregressive- 
moving average (ARIMA) model. The model is often written 
ARIMA(p,d,q), where p and q retain their previous meanings, 
and d represents the order of (non-seasonal) differencing. 
The model may also be generalized by including an additional 
constant term 84 that will have the effect of adding a 
deterministic polynomial trend of order d. The value of aS 


may be found in general by the formula 
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where 


The general ARIMA(p,d,q) model can be written compactly as 


d a 
(1) o(B) V z. = a + 9(B) a, 
An equivalent form for ARIMA(p,d,q) models that is 


sometimes seen is called undifferenced form: 
WKB) sz = “Gees O(B)i ca. ; 


where 


The "differenced" form, equation (1), is the form usually 
seen. In the differenced form, the transfer function $¢(B) 
is assumed to be stable; that is, all the roots of the 
difference equation $(B) = 0 are outside the unit circle. 
Giceariy, ¥(B) 1S not a transfer function of a stationary 
series when d > 0, since ¥(B) has dad roots on the unit 


eircle. 


RS, 


4. Seasonal ARIMA Models 

The general class of ARIMA models previously dis- 
cussed can be further generalized to allow for the modeling 
of seasonal models with the addition of appropriate seasonal 
operators and parameters. The seasonal model is called 
an ARIMA model of order (p,d,q) x (P,D,Q) ., where p,d, and 
q retain oie previous meanings, and P, D, Q and s are as 
defined in Section A of this chapter, page 17. The ARIMA 
(p,d,q) x (P,D,Q), model can be written as the product of 


nonseasonal and seasonal operators: 


(1-9,B) (1-¢) 158 VV 
Es THE AUTOCORRELATION AND PARTIAL AUTOCORRELATION FUNCTIONS 
Mathematically, the autocorrelation at lag k, denoted 


Pee is defined as: 


L El (z, - H) (2,4, Saat 
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The autocorrelation function of a time series describes the 


association (mutual dependence) among values of the same 
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series taken at different time periods, the difference in 
time being referred to as the lag, denoted k. Autocorrela- 
tion of itself implies nothing about a change in one 
variable causing a change in another. However, the auto- 
correlations may provide important information about the 
structure of a data set and its pattern. Ina set of com- 
pletely random data the theoretical autocorrelations among 
SuccesSive values will always be zero, whereas data values 
of strong seasonal or cyclical behavior will be highly 
autocorrelated. This P, can be estimated using the time- 


averaged sample autocorrelation, denoted r defined as: 


peu 
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A plot of the autocorrelation function versus the lag k, 
called a correlogram, is very useful for the purpose of 
determining if a process is stationary and for identifying 
the appropriate model. 

Another function useful in identifying the appropriate 
model for a given series is the partial autocorrelation 
function, often called the pauto function. Partial auto- 
correlations are analogous to autocorrelations in that they 
indicate the relationship of the values of a time series to 
various time-lagged values of the same series. However, 


they differ from autocorrelations in that they measure the 


oh 
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strength of the relationship between values of the series 
of various lags after the effects of other lags have been 
removed. In effect, they show the relative strength of the 
relationship that exists for varying time lags. For time 
series models of the types to be considered, the partial 
autocorrelation coefficients can be calculated several 
ways. One method of calculation, described by Box and 
Jenkins [Ref. 4, Sec. 3.2.5] begins eb the following 
equation, satisfied by the autocorrelation function, where 
the dy are autoregressive parameters: 
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This leads to what are known as the Yule-Walker equations, 


which may be written 
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In general, for Pre! the determinant in the numerator has 

the same elements as that in the denominator, but with the 
last column replaced with Pye The quantity dee an luncteion 
of the lag k, is called the partial autocorrelation function. 
Other methods for calculating partial autocorrelations include 
successively fitting autoregressive processes of orders 
mee;e---, K by least squares and picking off estimates of 
Oe a tee che of the last coefficient fitted at each 
stage. Another method, discussed in Appendix A3.2 of 

Ref. 4 and due to Durbin, generates estimates of an auto- 
‘regressive process of order k+l recursively from estimates 
of autoregressive processed orders k and less. It is 

derived by observing recursive relationships from the 


Yule-Walker equations. The recursive formulas are: 
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Another way of looking at partial autocorrelations 


is to consider a time series {z,} and the model 


Z = b. + b-2Z + DZ steve: SP ee ie ; 


where: the b values are the least squares estimates of the 


linear regression coefficients (8's) in the model 


Let ant be the residual of zy after removing the linear 


effect of z such that zZ i ey pe 


t-1! °°°! 7t-k+1 t’ t t t 
Then, the partial autocorrelation of lag k, denoted et, 


from Zz 


is defined to be the simple autocorrelation of lag k for 
the adjusted series (2,125, ese z,)° 

The exact expression for partial autocorrelations of 
a moving average process is complicated, but an approximate 


duality between autocorrelations of an autoregressive 
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process and pautos of a moving average process implies 
that the pautos of a moving average process would have the 
same general behavior as the autocorrelations of an auto- 
regressive process of the same order. This approximate 
duality also implies that the pautos of an autoregressive 
process of order p will have the same general behavior as 
the autocorrelations of a moving average process of order 
p. Thus, one function can be examined and used to confirm 
the other, and may be of use in model identification, discussed 
in section G of this chapter. 
G. COMMENTS ON MODEL SELECTION EMPLOYING AUTOCORRELATIONS 

AND PARTIAL AUTOCORRELATIONS 

The purpose of this section is to provide a brief over- 
view of the methodology for selecting the proper model for 
a given time series, through both examination or prior 
knowledge of the data itself and examination of plots of 
the autocorrelation and partial autocorrelation functions. 

A not yet stationary series should first be differenced 
or transformed until it appears to be stationary; character- 
istic of the autocorrelation plot (correlogram) of a non- 
Stationary series is a very slow damping out of the auto- 
correlations. When this property of the correlogram is 
observed, the analyst should difference the series until no 
further improvement is reached. For ARMA(p,q) models, the 
tentative identification of the model class, that is, the 


determination of p and q, can be done by examining the 
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sample autocorrelation and sample partial autocorrelation 
functions of the given time series with the theoretical auto 
and pauto functions of members of the general linear class. 
For most stationary time series, an adequate fit can be 
found in a model with p and q relatively small, say three 
or less. 

The following figure displays examples of autocorrela- 
tion and partial autocorrelation versus lag plots for several 


basic model types, described in Ref. 17. 


Figure 3. Examples of Auto and Pauto Plots 


Figure 3(a). Auto and Pauto vs. Lag for AR(1) Model 
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Figure 3(b). Auto and Pauto vs. Lag for AR(2) Model 
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Figure 3(c). Auto and Pauto vs. Lag for MA(1) Model 
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Figure 3(d). Auto and Pauto vs. Lag for MA(2) Model 


Auto Pauto 
1 as 
0 kK 0 [ k (time lag) 
ait -1 

Auto ol Pauto 
ih iil 
0 ce) 0 k (time lag) 
- -1 


38 


ctuet bee ost 


Figure 3(e). Example Auto and Pauto vs. hag FOr. 
Mixed ARMA(1,1) Model 
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As can be seen from the figures, there is a wide variety 
of possibilities for appearance of the correlogram and 
partial autocorrelation plots for even the simplest of 
models. Sometimes the pattern of estimated autos and pautos 
can be quickly classified in terms of one of the simpler 
basic models; however, most real data will generate auto 
and pauto plots that resemble those of Figure 3(e), those 
of a mixed model. Here some free association may be re- 
quired to infer a pattern from the correlogram, or more 
fehan one pattern may be implied. However, most of the time 
a selection of something is possible. As is evident, this 
is a somewhat subjective process, where the quality of 
identification will improve with experience. 

The following table, adapted from Box and Jenkins 
[Ref. 4], summarizes the properties of autoregressive, 
moving average and mixed ARMA processes; an understanding 
of its contents will be helpful as thumbrules for initial 
model identification. 

It can be noted from examination of Table II that there 
exists what can be termed a "duality" relationship between 
-the autocorrelations and partial autocorrelations of pure 
autoregressive and moving average processes. For example, 
Bae plot of autocorrelation versus lag for an AR(p) process 
would appear the same as the plot of partial autocorrela- 
Exom versus lag for an MA(q) process, for p = q. Similarly, 


the plot of partial autocorrelation versus lag for the AR(p) 
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Process 

Model in terms A 

of previous (B) 2, =a, 

Zz, Values 

e 

Autocorrelation Infinite and tails 

Function off; composed of 

damped exponentials 
and/or damped sine 
waves 

Partial finite; that is, 

Autocorrelation there will be p 

Function non-zero partial 
autocorrelations 

Table II. 


Autoregressive 


Moving Average 


Process 


finite; that 
is, there will 
be q non-zero 
autocorrelations 


Infinite and 
tails off; 
dominated by 
damped exponen- 
tials and/or 
Sine waves 


Mixed 
Process 


$(B)Z, = 6(B)z, 


Infinite and 
tails off; 
composed of 
damped exponen- 
tials and/or 
damped sine waves 
after the first 
q-p lags 


Infinite and 
tails off; 
composed of 
damped exponen- 
tials and/or 
Sine waves after 
the first p-q 
lags 


Properties of Basic ARMA Processes 


process would appear the same as the plot of autocorrelation 


versus lag for an MA(q) process, where p = q. 


For seasonal series, where differencing has been unable 


to remove all the nonstationary seasonal behavior, the auto- 


correlation and partial autocorrelation plots will generally 


exhibit a "Spike" at lags equal to integer multiples of the 


period of seasonality. 


For example, 


with a period of 12, then the autocorrelations 


if the data are seasonal 


Tio: ee and the partial autocorrelations 

12,12" 24,24! 36,36’ --.- would be amplified. This would 
indicate a need for one or more seasonal parameters to be 
included in the model. The following figure is an illus- 


tration of such an auto/pauto plot. 


Figure 4. Example Auto and Pauto versus lag for mixed 
ARMA(p,q), with seasonality remaining of 
period s = 6. 
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H. ESTIMATION OF MODEL PARAMETERS 


Suppose a time series has been tentatively identified 


as an ARIMA(p,d,q) model: 


estimated, (db) 4b5 Pep LOL pO 8 
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Box-Jenkins procedure separates this estimation into two 
parts. First, estimates are obtained for the AR and MA 
parameters $= (d41%5" has d= )imandoe =n (se 7 yamine 
then estimates are calculated for 8G and oe which are func- 
tions of $ and 6. 

The usual procedure is to select those parameter values 
$ and 6 that minimize the sum of Squared model errors 
(residuals). Let ee af - 4 - 5 AR ie e,,) and 

Pa . 
W, = V (Z, Ws Then, it can be shown that Peo E(w] 


and that the model can be rewritten as 


or 


ae gee (8) 0B) (W, = 19 


Now, let W, the sample mean, be the estimate of Le (i= 


a> 0, then oe is wsually zero): Also, Let 
Sess (EOE) (We W) 
e t ; 
Set 
N 
Sto ,eva ee as 
ee t ! 
t=1 


where N = n-d is the length of the differenced series. The 
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objective is to select those parameters $ and 8 such that 
S($,8) is minimized over all values of the (p+q)-dimensioned 
parameter space. Since the equation in S is nonlinear 

in the parameters, iterative search methods are generally 
used to determine optimal model parameters. Having found 
these optimal parameters and the minimum S, of can be 
estimated using 


ae = $§($,6)/N-p-q . 


Chapter III contains a more detailed description of a 
nonlinear algorithm to estimate parameters of a generalized 


seasonal or nonseasonal ARIMA model. 


I. DIAGNOSTIC CHECKING OF MODELS 

After the model has been tentatively identified and 
parameter estimates have been calculated, the next task is 
to test whether or not the original model specification 
was correct and that the model itself is adequate in fore- 
casting power. The process of testing the model can take 
many forms, but will usually include at least the following 
two steps: 

1. Generate a simulated series from the estimated 
model and compare the simulated series and its 
auto and pauto functions with the original series 
and its respective auto and pauto functions. This 


comparison is essentially subjective. 
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2. Calculate the residuals of the estimated model, the 
a.'s, and compare the properties of the residuals 
with those assumed for the shock terms of the actual 
process. The residuals should be normally eset 
buted and uncorrelated with each other; that is, there 
should be no discernable structure in the residuals. 
There are many quantitative statistical tests and 
data analytic tools that can Be tapemnea to the 
residuals to test hypotheses of normality and zero 
autocorrelation. 

A plot of the autocorrelation and partial autocorrelation 
functions of the residuals will provide not only a test of 
whether or not the residuals are uncorrelated, but, if they 
eee Correlated, the plots can be used to suggest improve- 
ments to the model. For example, suppose the model was 


tentatively specified as the ARMA(1,1) model below: 


(R= LOTS) ize oa 22) SWS CL se 20. 7B) ay 


Suppose also that the autos and pautos of the model residuals 
= u,_, where the u_'s 


c € c 


: sh ; ; 2 
are white noise, uncorrelated with variance oe. Then, 


Suggested the model (1 - 0.3B) a 


these two models can be combined into the ARMA(2,1) model: 


pia Cass) a (ies s5B) (ze —e2) = 20 + 0.7B) a 


te ec 


This should be a refinement over the original model. 
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J. FORECASTING 
The objective in forecasting is to predict future 

values of the time series with as little error as possible. 
The criterion most often used for forecasting is to calcu- 
late that forecast which minimizes the expected mean square 
forecast error. Therefore, if 2, (2) denotes the forecast 
for lead time 2% from time origin t of the value Ze ag! 
the objective is to find Zz, (2) such that the objective 


function 


is a minimum. This forecast is given by taking Z, (2) as 


given z Z 


the conditional expectation of z Looccer 24? 


t+2’ 


2 02) Caan Zz 


Le rl cig ng veer Zz] 
The forecast can be easily generated recursively from the 


Mathematical model, utilizing the facts that 


ElZ,_4|Zp7+++ +2] Rees 
Bet 7 = 0,1,2... (where t is the current time), and 
BZ 4lZpr---72y] = z,(j), BOTT lg eigk xh ase 


46 


ove 
ae . 7 a 
' 


: ‘ 
a 
sotbesy) oF af praaeas 


so1re alsttl’ 26° eae 


are! ‘_s 722 ,s rote zene 


oo : ~ 3 rf 12 a ie adie + moss ’ 


at (t+ of wl 
Live i ’ oes = 


and 


0 for times > t 
Ela, ] 


a forptimes..< t 
For example, consider the Box-Jenkins model 


ia = (0: sBat 0.6B°) Z =— (1 + 0.3B) a 


€ € 


which could be written in the expanded form as 


Now, assume it is known that Z100 rSo.L44, Z99 isjls0e and 
the calculation has been made to obtain the residual 

A100 = 0.2. Then the forecasts of Zi0L" 2402 and 2403 
made from an origin of t = 100 can be calculated from the 


model as follows: 


Linse, tor £ = 1, set 


ElZ1 9112199" 2g97°° 124) | 


Cy 9S El0.52, 9 = 0.624, + Ai01 + 0.38, 9) 


’ 


and 
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zeae SphOb7 Seor6o4u00.06 


Il 
fo) 
e 
bh 
fon) 


Next; for 2 = 2, set 


Z199(2) = ElZygol2i 91/21 00/299" 24) 

agile) °= E[0.5Z, 9, SuOeGz f S495 * ©: 28400 
ines? 711 ; 

ee OLS) (2 (1) = (0.6) (2544) + 0 + 0 4 

Zale) = Os) (Oct6) = (0.6) (lia) = <0276. ; 


Scimtlariy, for % = 3, set 


21003) = El219312102°%101"7100,299' "21! 
Zareots) = E[0.52Z, 45 - 0.62, 5) a aie Oday 

123 99" 2997°* 712! 
gos) = (GES) (z4 4 (2) (Woqvie ahh + O40, 
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and then 
(Ss) =) (075) (=—0-2 76) = (0.6) (0-16) = -0.48 . 


Clearly, this process can be continued into the future as 
long as desired, recognizing that the expected forecast 
accuracy will decrease as 2 increases. 


Let 


denote the forecast error for lead time 2 beyond the forecast 


Eeime Otigin t. It can be shown that e,_(2%) is given by 


t 


the relation 


+ 


* Mpnateer ’ 


SD Feo Gea PY Pee pa 


where the scalar weights Ve are determined from the equation 
i =d 


Baas woe Aiba = B)iy. 6 (B) 


using the known parameter values. The variance of the 


forecast error is given by 


2 - 2 2 De Did 
fee) = (Ly) tw. + Wj) (ol) = {1 + ) eae 
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From these relationships, a confidence interval of n 
standard deviations about a forecast of lead time 2% would 


be given by: 


c— 1/2 


A = “Aa hk 2 A 
eezett))) = -2.(2) 2 n[1 + } v5) ie 


zie 


The value of n could also be though of as the percentage 
point for the desired level of confidence using the Normal 
O@estrabution [Ref. 4]. It can be noted from the definition 
of e, (2) above that the one-step-ahead forecast error, 


e, (1) is simply the value a that 1s, 


+1’ 


This is an explanation for the use of the term "residual" 
to refer to the white noise or shock terms. Additionally, 
from the foregoing it is evident that the forecast error 
variance is a non-decreasing function of the length of 
the forecast lead time 2; therefore, the confidence bands 


must become wider as the forecast lead time increases. 
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III. GENERALIZED ARIMA MODEL PARAMETER ESTIMATION 


The purpose of this chapter is to describe the compu- 
tational details of the algorithm WMAROQRDT, employed with 
appropriate subroutines to estimate Box-Jenkins parameters 
of a generalized seasonal (or non-seasonal) ARIMA model. 


This generalized model can be written parsimoniously as 


follows: 
s - s 
ONE) (ES Oar, ST THOS TOMB POKB™ ya, 11 
where 
: aD... 
W, is V Vazer 
and 


zy is the transformed time series (for example, 


Ve oS 
Zz in] 2 or Z, vZ,) - 


The model employs a Marquardt-type [see Ref. 13] non- 
linear least squares algorithm for determination of the 
model parameters, Searching in the parameter space for the 
set of parameters which minimizes the sum of the squared model 


residuals. 


5 
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A. MATHEMATICAL DISCUSSION OF THE MARQUARDT ALGORITHM 
Reference 13 describes Marquardt's methodology for 
development of a general algorithm for ee esti- 
mation of nonlinear parameters. Most algorithms for least- 

Squares estimation of nonlinear parameters have been cen- 
tered about either a pure linear iterative model based on 

a Taylor series expansion or using some form of the method 
of steepest descent (maximum negative gradient). Since 
both methods have severe potential pitfalls (Taylor series 
due to divergence of the iterates due to an unfortunate 
choice of initial values, the steepest descent due to slow 
convergence after the first few iterations), a method termed 
"Maximum neighborhood" method was developed by Marquardt 
which, in effect, performs an optimum interpolation between 
the Taylor series method and the method of steepest descent, 
the interpolation being based upon the maximum neighborhood 
in which the truncated Taylor series gives an adequate 
representation of the nonlinear model. 

As discussed in Ref. 14, the idea of Marquardt's method 
can be briefly explained as follows. Assume the algorithm 
begins with a vector of initial parameters 0G GO ee tea) 
If the method of steepest descent is applied, a certain 
vector direction nee where the subscript g represents the 
gradient (which gives the direction in which the rate of 
change is the greatest), is obtained for movement away from 
the initial point values in the parameter space. Due to 


potential nonlinearities in the sum-of-squared-residuals 
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surface, S (5795795794) » in the parameter space, hs may 

be the best local direction in which to move toward opti- 
Reity (minimal sum of squared model residuals), but may 
not be the best direction. However, the best direction 
must be within 90 degrees [Ref. 8] of hos or else the 

value of S(o,%,8,9) will increase locally. Now, the 
linearization, or Taylor series method, may result ina 
different correction vector, h,- Marquardt found experi- 
mentally that for a number of practical problems, the angle, 
say &, between the vectors he and hy, fell in the range 

80 degrees < € < 90 degrees. In other words, the two 
directions were nearly normal in most cases. The Marquardt 
algorithm then, provides a method for interpolation between 
the vectors Ae and hy employing a sutiable step size, which, 
in general, converges more quickly than either of the two 
"parent" methods alone. The remainder of this section des- 
cribes the mechanics of the algorithm, as applied to the 


estimation of Box-Jenkins parameters. 


B. MODEL INPUT INFORMATION 


The following is the required input information to 


the WMAROQRDT parameter estimation model: 


pws the transformed time series (for example, 
Qo es ‘= 
Wy ln Zi Or WwW, v2, 
N the length of the time series 
s the length of the seasonal period (for no 


seasonality, s = 1) 
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eerie. 


~f3c0 bene’ syom OF . ni’a 
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——<_ + 1 + 1. | 


P the number of seasonal AR parameters, P > 0 


p the number of non-seasonal AR parameters, 
p> 0 
Q the number of seasonal MA parameters, Q > 0 
g the number of non-seasonal MA parameters, 
q> 0 
NDIFNS the number of non-Seasonal differences to 
be taken 
NDIFS the number of seasonal differences to be taken 
6 = (® Pieralcrg ) initial estimates of seasonal AR 
Ls P,s 
parameters 
$ (by rbgr---0d ) initial estimates of non-seasonal 
BP AR parameters 
6 = (6 ,0 abet hte ) initial estimates of seasonal 
NS Q,s MA parameters 
6 = (6,85, i 84) initial estimate of non-seasonal 


MA parameters 


C. ALGORITHM DESCRIPTION 

Initially, the model accepts inputs of the time series 
length and the (transformed, but undifferenced) time series 
itself from the user's disc, where it is stored in File 
FTO2FOO1. The remainder of the input data is then entered 
interactively on the terminal as prompted by the model. The 
model then forms the initial parameter estimates into vectors 
St lengths p, P, q, and Q, and employs them to calculate 


an initial model residual sum of squares. 
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Prior to initiating the sum of squares calculation pro- 
cess, the model unravels the polynomial of AR and MA operators, 
and forms a vector of coefficients for both the AR and MA 
sides of the model. For example, consider the relatively 


simple (2,1) x (2,1)55 ARIMA model, written 


2 a2 ie a - 2 cc 12 
(1-9, B-$.B )(l - 1128 Wy Sy tonlJli 6,8 6B (2 ©, 128 )a,- 
Unraveling, the model becomes first: 
2 or Ley seo . 
ao B-$58 )(W.-o) 15,35) = (1-0,B-05B°) (a, -©) 1244-12)" 
and then takes the form: 
ees 1 28 i ear Pa po e137 at es27 22112" e=04 
ee Pea pt Pee tea) 3k’ 2 ea? 2°, 12"e-14 . 
This is represented in the computer as four vectors: two 


vectors, 1 and y, composed of the expansion AR and MA 
coefficients respectively, and two vectors composed of the 
indices of the time series terms associated with the coeffi- 
cients. For the model in the example above, these vectors 


would be written as follows: 


T = (1,9) gor dyr (Oy MR 79) bor (b9) () 42) 
AIDE peek neh 2um ks pL31.2,414) 
Sma oe anio, tere (Cy, ap) ear (SKE 49)? 
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This format for unraveling of a generalized product of 
polynomials of seasonal and non-seasonal parameters permits 
simplified calculations in the determination of the sum of 
squared residuals, since the arithmetic need only be per- 
formed for the non-zero parameter values, those contained 
in the 7 and y vectors. In the actual calculations, per- 
formed in FORTRAN, the index vector origin is translated 
to the right by the amount p + Ps + q + Qs, called IADDIT 
in the model, in order to prevent obtaining zero or negative 
subscripts. 

The calculation of the residual sum of squares is a 
three stage process. First, the model calculates ten 
values of e,, the "forward" white noise terms, using the 


c 


following relationship: 


The assumption is made that (e 7+ 72y) are all equal 


Peacno- 


to zero. Having found the desired values for (C19 7egr ee ey) 


we can now backcast the necessary values of W in order 


ter” 
to enable us to calculate the required white noise terms 


a as follows: 


t! 
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solving for (W\,W_,/W_5,---/W ). Now, with the 


1- (p+sP+q+sQ) 


W.'s backcast, we can proceed to calculate the estimated 


white noise terms (residuals) using this relationship: 


The values of a, are’ Calculated for t = 1, .«.., Ns. Once the 


a's are calculated, it is a simple matter to calculate the 


sum of squares, using 


2 


N 
S(o,6,9,9) = ) ee ° 


t=1 


Having calculated the residual sum of squares for the 
initial parameter estimates, we form the parameters into 
a vector 8, such that 


8 = (8,785, ceere By) ’ 


where k = p+P+q+Q ; that is, 


Syd) 


(Qeep-de- 


FEL UGLED 


The subroutine PARSH then calculates the derivatives 


Ja 


if ea oat 
Lt Of a. 
a 
Syer all values of t = 1,...,N and i=-l1,...,k. Using the 


residuals {a, } calculated earlier, the derivatives are 
estimated numerically using a perturbation of each parameter 


of amount 6, such that 
a, — fa, (By r-- +r By) ay (By ree +1 Bit5,---18)) t/85 


Then with {a,} and x5 Supplied for the current parameter 


rt 


values, the following quantities are formed: 


meethe kKxk matrix A = {Aj 43 , where 
N 
ni4 7 ” Ce ae 
=1 
weeeethe vector G = {g;} , where 
N 
Gg. °= 4 (x, .) (a,) 
t=1 
Se the scaling quantities Ds = VAS . 
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Then, the modified (scaled and constrained) linearized 


equations 


LfAsimeine} oS tig¢g*] 
are formed, where 
a = Bae Bees yy 2 F j A 
A;,* = A;,* emuee ge nee! OT. oe 
and 
Gi, Ves ,98/ Bin. 


It can be noted here that the beret! Matrigeis wcofact ca 
correlation matrix for the model parameters. Those equa- 
tions are then solved for h*, which is scaled back to give 


the parameter correction vector Sl where 


Renee ee /D 
a 3 / J 


Then, the new parameter values are formed from the h vector, 


such that 


a) 


and the sum of squared residuals for the BNEW set of 
parameters is calculated. 


Now, if S (By oy) Soe ), the parameter corrections h 


OLD 
are tested. If all elements of the h vector are smaller 
than some epsilon, say .0001, then convergence is assumed 
and the kxk matrix [ax] is used to calculate the covari- 
ance matrix of the estimates; otherwise, Borp is reset to 


BNEW’ the constant € is reduced by a factor F, say .l, and 
computation returns to the calculation of a new set of 
derivatives. 

However, if S (Bey) > S(Borp): the constraint parameter 
€ is increased by the factor F, and computation is 
returned to the stage where the A matrix is formed. In 
all but unusual cases, a minimal sum of squares will be 
found. However, an upper bound can be placed on &, and 
when it is exceeded, the search is terminated. When con- 


vergence is reached, the desired output information is 


calculated. 


D. MODEL OUTPUT AND DIAGNOSTIC CHECKS 


The output includes: 


H p= (Byr-++7 By) the least-squares estimates of 
O7r- <1, 
Presse Py 
Sy ea 
Oyre- 195 
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Zs oS the residual variance, N-p-q-P-0O S(@), 0, o1,0) 
sz. V the covariance matrix of the parameter 
estimates, formed from 
rf os =i x2 
V= eee = (AA) | ce 
4. Si the standard errors of the parameter 
estimates, 
Si = UES 7 1 = SC HiOare ag sel Ole 
5. Ri; the correlation matrix, obtained from 
Reka" VN 
25 cs 5 gp ge 
6. 6, the overall constant term, where 


n~ 


oe =u G , £or 


E> 
ll 
Ale 
= 


Finally, diagnostic checks on the model are performed. 
This includes calculation of residual autocorrelations 


ro ik), obtained from these formulas: 


rar(k) = ea tie eaA(0) - 
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where 


Be ic cas 
eal 


and k goes from 1 to the maximum desired lag K, usually 
40. Additionally, a chi-square statistic is calculated 


from 


K 

Z a 

= N yen i) F 
k=1 


2 
X(y) 


and is compared with a chi-square distribution with 
m— rk —- Pp- gq- P-Q degrees of freedom. 

Examples of input methodology and appearance of output 
for this model are contained in later chapters of this 


paper. 
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TV. DESCRIPTION OF THE TIME SERIES EDITOR 


In this chapter, descriptions are given for each of 
those programs in the Time Series Editor which interact 
with the user. There are fifteen separate program modules 
currently included in the Editor; there are other programs 
utilized by the Editor that are completely transparent 
to the user, and are not described here. However, listings 
for all programs contained in the Editor are provided in 
Appendix D to this report. 

No attempt is made in this chapter to describe the 
actual mathematical calculations or algorithms that are 
performed by the programs. Rather, the objective is to 
provide general descriptions of what each program can do 
for the user, and how the user interacts with the programs. 
The programs described in this chapter are, in the order 
presented: TIMESER EXEC, ZFORMAT, CMSWORK, TRANS, DIFF, 
PLOT, AUTO, ESTIMATE, YESTSEAS, WMARQRDT, XSUMSQ, FORECAST, 


ROOTS, SIMULATE and GENERATE, and HELP. 


A. THE EXECUTIVE PROGRAM 

The heart of the Time Series Editor is a master program 
@artea TIMESER EXEC which provides file control for all of 
the other program modules, controls input and output, super- 
vises the necessary CP/CMS protocol and provides instructions 


interactively to the user concerning the contents of the 
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Time Series Editor and how each program can be used. 
The TIMESER EXEC program is written in the IBM-360 CP/CMS 
Executive language. It is the only program in the package 
not written in FORTRAN; consequently, it should be the only 
program that would need modification if the Editor were to 
be adapted to another FORTRAN-capable time-sharing system. 
After the user has logged into CP/CMS and linked to 
the disc space containing the Time Series Editor (instruc- 
tions for this procedure are provided in detail in the 
User's Guide to the Time Series Editor, included as Appendix 
A to this report), the entire Editor package is made avail- 
able to the user by the command TIMESER. On entry of this 
comand, the executive routine, TIMESER EXEC, can provide 
a guided tour through the Editor. It will briefly des- 
cribe what the Editor can provide and asks the user what 
Option he wants to use. On the basis of the user's response, 
the Editor then advises the user what input data is required 
and how it should be entered, either through a data file 
entered offline or parametric data entered interactively 
at the terminal. When the user selects an option for exe- 
cution, the TIMESER EXEC routine loads the appropriate pro- 
gram package and automatically manipulates the required 
input and output data files. For the more experienced user 
who does not require detailed user instructions, there is 
a shortened version of the TIMESER EXEC. Upon entry into 


the shortened version, the user is immediately asked which 
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option he desires. This version of the Editor is entered 
by adding the argument "s" when logging into the Editor; 
that is, type TIMESER S when logging in. The TIMESER S 
(shortened) version of the Editor provides exactly the 
same program options to the user as the longer TIMESER 


version. 


B. DATA AND FILE MANAGEMENT 

Whenever data is required, the user is prompted at 
the terminal by either the TIMESER EXEC or the program 
module being executed. In most cases, the necessary user 
response is a short alphanumeric character input during 
execution using the terminal keyboard. However, when the 
time series data itself is entered, the user may enter the 
eaeariy Card deck offline. Similarly, most output is pro- 
vided right at the user's terminal. Output such as listings 
of transformed series or plots are sent to the offline 
printer to conserve time and provide hard copies of the 
results. Detailed descriptions of the input and output 
requirements of each program module are provided in this 
chapter, and also in abbreviated tabular form in Appendix 
A, the User's Guide. However, there are some general prin- 
ciples of data management that apply to all programs. 
These are described in this section. 

me OLt line Data Entry 

The Time Series Editor requires that the user's 


time reside in FILE FTO2F001. The data can be read offline 
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via cards and directed to FILE FT02F001. An example of the 
required format forthe input card deck for the Time Series 
Editor is provided in Figure 5. It is noted here that the 
data format expected by all programs is FORMAT(5F15.6), 
preceded by a card containing the length of the time series 
in FORMAT (I3). 
2. ZFORMAT Program 

In many cases, the user will want to analyze time 
series data that are already available ina format other 
than that required by the Editor. To spare the user the 
tedious task of retyping his deck in the required format, 
the ZFORMAT. program was written. seg program converts a 
time series data file (FILE FT03001) written in any FORTRAN 
eeemat Into the proper FORMAT(5F15.6) onto FILE FTO2FO001. 
The length of the series is entered as an input on the 
terminal. The data are then ready to use in the Editor, 
waeethne Original file left in FILE FTO3F001. Figure 6 
shows a sample input deck for such non-standard data. 


Figure 7 gives a sample session where ZFORMAT is used. 
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timeser s 

ENTER LETTER FOR OPTION YOU WANT. 
ts YOUR DATA IN FILE FTO3F001? 
EXECUTION BEGINS... « 


iS YOUR TIME SERIES DATA NOW IN FILE FTO3F001? 
¥ 


ENTER THE LENGTH OF YOUR TIME SERIES VIA FORMAT I3. 
010 


NOW ENTER THE FORTRAN FORMAT FOR YOUR TIME SERIES DATA, 
ENCLUDING PARENTHESES; FOR EXAMPLE, TYPE: (5F15.7). 
pO 2 


TOUR FORMAT IS: 
FORMATF10.2 

fi LHS CORRECT? 
n 


NOW ENTER THE FORTRAN FORMAT FOR YOUR TIME SERIES DATA, 
PNCLUDING PARENTHESES; FOR EXAMPLE, TYPE: (5F15.7). 
feeer 0’. 2) 


YOUR FORMAT IS: 
FORMAT (F10.2) 

fo LHIS CORRECT? 
y 


YOUR DATA IS NOW TRANSFORMED INTO THE PROPER FORMAT 
Bem Usk IN THE TIME SERIES EDITOR, LOCATED IN FILE 
BPrO2FO001. YOUR ORIGINAL DATA IS IN FILE FTO3FOO01. 


DO YOU WANT TO GO AGAIN? 
n 


CONTROL RETURNED TO CMS 
R; 


Figure 7. Sample user session with ZFORMAT program 
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3. Interactive Keyboard Input 


a. Numeric Input 

When the user is prompted at the terminal to 
enter such numerical values as the number of AR parameters, 
estimates of parameters or the number of differences to be 
taken, he must enter these values before program execution 
can continue. The user should enter the data according to 
standard FORTRAN practice. For example, integer data should 
be entered (without a decimal point) for counts, Such as 
the length of a seasonal period, and for names beginning 
with the letters I through N; floating-point data, such as 
initial estimates of parameters, should be entered with a 
decimal point. Because typed decimal point overrides a 
floating-point format, it is not necessary for the user 
to concern himself with the format for floating-point input 
data. However, care must be taken when entering integer 
Gata because it must be right-justified in its format field. 
The user is advised in each case when the format is other 
than Il. For example, suppose the user desires to perform 
analysis on a series having only 90 observations. If the 
program that he is executing requires the length of the time 


series, the user will receive the following request: 
ENTER LENGTH OF THE TIME SERIES VIA FORMAT I3 


The user would then enter: 
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column 123 
b90 (where b represents a blank space). 


If the blank were omitted, the program would read the 
series length as 900, and problems in execution would 
mpocur. 
b. Alphabetic Input 

The Editor often asks the user to respond with 
alphabetic input. Generally, this is in the form of the 
response to a question requiring a yes or no answer, or a 
Bee eonse to the question of which program to execute. The 
editor has been programmed to read only the first letter 
of such responses. Therefore, the user need only enter the 
first letter of the response for each such inquiry. For 
example, the’ user should enter Y for yes, N for no, P for 
the PLOT option, X for the XSUMSQ option, W for the 
WMAROQRDT option, etc. The other alphabetic input is 
generally entry of titles to plots; these can be any 
combination of numeric and alphabetic characters, as long 
as they do not exceed 72 columns in length. For example, a 


Suitable title would be: 


PLOT OF AUTOS AND PAUTOS FOR SERIES C DATA (2 ORDINARY DIFFERENCES) . 


4. Output 


Most of the results are written out at the user's 


terminal. However, in some cases the output is printed 


fab, 


inetd & 
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offline to conserve time. Such results as plots and 
transformed time series are written onto various files 
PeTOSFOO], FTOSFO01, FTO8FO01L), and the TIMESER EXEC pro- 
gram causes them to be automatically printed offline under 
the user's identification number (USERID). The plots and 
files can also be printed out at the terminal using the 
usual CMS commands in the CMS mode, outside the Editor, or 
in the CMSWORK program (see next Aa However, the 
user should be aware that plotting graphs at the terminal 
is generally mutes a slow process. 
5. CMSWORK Program 

This program permits the user to enter the CMS 
environment to perform routine CMS functions without leaving 
the Time Series Editor. Commonly used CMS commands include 
file name alteration, file erasure, file offline printing 
Bespunching, obtaining a disc status, and listing all disc- 
resident files. However, essentially all CMS commands can 
be issued while in CMSWORK. Figure 8 provides an example 


of the use of the CMSWORK program. 


C. TRANS PROGRAM 

The TRANS program takes a given time series in File 
F102F001 and performs a transformation upon each data ele- 
ment. The options include a change of scale, a power trans- 
Permation, a shift of the origin, and a natural logarithm 


transformation. The original data is written onto a file 
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ENTER LETTER FOR OPTION YOU WANT. 
Cc 
ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 
stat 
fees) 2OS0eFILES; 263,REC IN USE, 33 LEFT (of 296), 89% FULL (2 CYL) 
fester * £tO2ZF002 
FILE NOT FOUND 
mee 15 LISTE * FTO2ZFO02 
Pitee{O0002) !'! 
sce * £tO02F001 
FILENAME FILETYPE MODE NO.REC. DATE 


SERG ELOZEFOO]L PL 3 97 16 
SERC ETO2ZE OOls Oeil 5 9/16 
LNSERG FTO2FOO1l Pl 3 9/16 
FILE FLOZEOOL, Pi 3 9/16 
ep gq £ 


Beebo. - NO RDR, NO PRT, NO PUN 

erase file f£t08f001 

Setlaine print file £t02£001 

Beeline print serg f£t02f001 

Sretine print serc £t02f001 

erase file ft02f001 

Stat 

feo): 26 FILES; 240 REC IN USE, 56 LEFT (of 296), 81% FULL (2 CYL) 
“&goto -ques 

DO YOU WANT TO GO AGAIN? 


n 
CONTROL RETURNED TO CMS 
R; 


Figure 8. Sample user session with CMSWORK program 
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called DATA FTO2F001, while the transformed series is 
written onto FILE FTO2F001. The parameters of the trans- 
form itself are written onto FILE FTO7FO01. The TRANS 


program is self-contained, using no other programs. 


D. DIFF PROGRAM 

The DIFF program takes a given time series in FILE 
FTO2FO001, and performs non-seasonal and/or seasonal differ- 
encing operations on the series in order to help achieve 
stationarity for proper modeling. The user must indicate 
if seasonal differences are required; if so, the length of 
the season is input. The user is then required to enter 
the number of seasonal and nonseasonal differences. After 
execution of DIFF, the original data are in FILE FTO2F001, 
and the differenced data are written onto FILE FTO3F001. 
To use the differenced data in any other program, CMSWORK 
program must be entered and the files altered to save the 
Original data (if desired) and to transfer the differenced 
data to FILE FTO2F001. The CMS commands to perform this 


shift of files would look like: 


ALTER FILE FTO2F001 Pl. SAVEDATA FTO2F001 Pl 
R; 

ALTER FILE FTO3F001 Pl FILE FTO2F001 Pl 

Ry 


The DIFF program uses the IMSL subroutine called FTRDIF. 
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E. PLOT PROGRAM ; 
The PLOT program plots any given time series which 
resides in FILE FTO2F001. Other than the time series, 
which must be on the user's disc or entered offline, the 
program requires only that an identification title for the 
plot itself be entered during execution. The plot is auto- 
Matically printed offline. The PLOT program uses a modified 
PLOTP subroutine called PLOT8 from the IBM Scientific 


Subroutine Package Library (SSPLIB). 


F. AUTO PROGRAM 

The AUTO program takes a given time series in FILE 
FTO2F001 and calculates summary statistics of utility in 
Box-Jenkins modeling. The statistics include the sample 
mean, variance, autocorrelations and partial autocorrela- 
tions for lags one though 40. Part of those statistics 
are printed out at the user's terminal, while plots of the 
autocorrelations and partial autocorrelations versus lag 
are printed offline. These plots will also be written onto 
the user's disc in FILE FTO8F001, and may be plotted at 
the terminal in CMS outside the Editor or in the CMSWORK 
program; recall that thisis a time-consuming process. The 
AUTO program provides the user with essential information 
about stationarity, seasonality and model identification. 
With the Box-Jenkins procedures, the second moments (auto- 
correlations and partial autocorrelations) are valuable 


tools in tentative model identification. Additionally, AUTO 
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is useful in the diagnostic checkout phase of model 
building, when examining the model residuals for structure. 


The AUTO program uses the IMSL subroutine FTAUTO. 


G. ESTIMATE PROGRAM 

After the user has tentatively identified a model for 
a series through analysis of the series plot, autocorrela- 
tions and partial autocorrelations, ESTIMATE should be 
executed to calculate maximum likelihood estimates of the 
model parameters. ESTIMATE is written to calculate parameters 
for non-seasonal models only; for seasonal models, the 
programs YESTSEAS and WMARQRDT should be used for parameter 
estimation. The ESTIMATE program requires that the series 
to be modeled (already transformed and differenced to 
achieve stationarity) reside in FILE FTO2F001; additionally, 
the program will instruct the user to input the number of 
MA and AR parameters in the model, as well as the number of 
differences taken. The general model for which ESTIMATE 


calculates parameters is: 


$43) (Wighsai@ bt 0(Blica 


where we = vont and Zi is the transformed value of Zy- 
ESTIMATE calculates estimates of the AR parameters for the 
undifferenced form of the model, the estimated MA parameters 
and MA constant term, the residual variance, autocorrelations 


and partial autocorrelations of the residuals, plots of 


76 


ae 7 nx: 


fshom to seedq so0spS 
arutouxjge 2ot eisphtdad tebom ame 
OTUATY ontscoudie. Tenz ons 


—y 
he TA ssitiscebt visvisedasd 26d see 


2) rol 25i%36 aft Io ateyiseaa ¢ 


oi 


s 
2“o.cseLlertocesce 4B 


- nixem stele 


ITAMITes See 


2Leet offs. \fee 
‘ee 
si sist t® enoive Seas 


~ 
A 


these autocorrelations and partial autocorrelations, and 

a Chi-square goodness-of-fit statistic for the model. The 
original series to be modeled remains in FILE FTO2F001 

and the residuals of the model are stored in FILE FT02F001. 
ESTIMATE uses the IMSL subroutine FTMAXL (a nonlinear 


gradient search algorithm) to estimate the model parameters. 


H. YESTSEAS PROGRAM 

The YESTSEAS program is used to calculate initial esti- 
mates of autoregressive and moving average parameters for 
a generalized seasonal Box-Jenkins model. These initial 
estimates are used as inputs to the program WMARQRDT, which 
then refines them to determine the parameter estimates that 
minimize the sum of the squared residuals. The general 


model assumed in YESTSEAS is 
9(B) o(BS) 7ty2z, = 9(B)o(BS)a 
s*t ra 


The YESTSEAS program requires that the seasonal time series 
for which initial parameter estimates are required is located 
in FILE FTO2F001; it should be transformed, if desired, but 
generally not differenced, since YESTSEAS performs the 
required (seasonal and/or non-seasonal) differencing. 
Terminal inputs required by the program include the number 

of sdasonal and non-seasonal differences to be taken (if 

data is already differenced, enter zero for these values), 


the length of the seasonal period and the number of each 
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type of parameter to be estimated (that is, the number of 
each of the seasonal and non-seasonal AR and MA parameters). 
The output of the model, initial estimates for the requested 
parameters, is printed out at the terminal. YESTSEAS uses 
the IMSL subroutines FTRDIF, FTMAXL, FTAUTO and FTARPS, 


as well as a modified version of FTMAPS. 


I. WMARQRDT PROGRAM 

The WMARQRDT program uses the initial parameter esti- 
mates calculated by YESTSEAS or other means as Starting points 
for calculation of parameters of a generalized seasonal 
Box-Jenkins model which minimize the sum of the squared 
residuals. The general seasonal model assumed by WMARQRDT 
can be written: 

$(B)o(BS)V°Voz, = 8) + 6(B)O(B°)a,. 

Details of the Marquardt-type non-linear least squares 
algorithm employed in this program are provided in Chapter 
III of this report. WMARQRDT requires that the series to be 
modeled reside in FILE FTO2F001; it should be transformed, 
if desired, but generally not yet differenced, since the 
program performs differencing. Inputs at the terminal 
include number of seasonal and non-seasonal differences, 
length of the seasonal period, number and type of parameters 
desired, and initial estimates of those parameters. The 


program then calculates the least-squares parameter estimates, 
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the standard error of these estimates, the MA constant term, 
the residual variance, the autocorrelations and partial 
autocorrelations of the residuals, the final model sum 

of squared residuals and a Chi-square goodness-of-fit sta- 
tistic (the portmanteau test). Additionally, it plots 
offline the autocorrelations and partial autocorrelations of 
the residuals. WMARQRDT uses the IMSL subroutines FTRDIF, 
FTAUTO, LINVIF and MDCDFI, the SSPLIB routine PLOT8, and 

the Time Series Editor resident subroutines PARSH, MARQRT, 
SUMSQ, SWAPB and FORMB. The calculations are performed in 


double-precision arithmetic. 


J. XSUMSQ PROGRAM 

The XSUMSQ program accpets an already transformed and 
differenced series in FILE FTO2F001 and general seasonal 
and non-seaonal Box-Jenkins parameter values for a general 
seasonal model from the terminal. The program then calcu- 
lates the residual sum of squares for these parameters, 
giving the user a feel for the "goodness-of-fit" the 
parameter estimates provide in modeling the given time 
series. XSUMSQ uses the Time Series Editor resident sub- 
program XSUMSQ. All output is printed at the terminal. 


The calculations are performed in double-precision arithmetic. 


K. FORECAST PROGRAM 
The FORECAST program uses the estimated Box-Jenkins 


(Seasonal or non-Seasonal) model to compute forecasts of 
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the transformed and differenced time series. FORECAST re- 
guires that the time series to be forecasted reside in FILE 
FTO2FO01 in transformed form. During execution the input 
required at the terminal includes the forecast origin, the 
numbers and estimated values of the AR and MA parameters 

of the model, the overall MA constant, the plot origin 
index, the maximum forecast lead time, the order of differ- 
encing in the model, and a level of SiGmiticance for the 
forecast confidence limits. The program needs the transfor- 
Mation parameters from FILE FTO7F001, created when using 
program TRANS, and performs the inverse transform to return 
Bae Gata to its original form. FORECAST also undiffer- 
ences differenced data. The program output includes the 
forecasts up to the given maximum lead time, the deviations 
from each forecast for the confidence limits, and plots of 
forecasts and confidence limits plotted offline. The 
FORECAST program uses IMSL subroutine FTCAST and the special IBM 
SSPLIB Subroutine UTPLT8. The plot itself and output data 
are written onto FILE FTO8FO001 on the user's disc; since 
this plot file takes considerable space, it is usually 

best to erase it after execution of the program is completed 


and the plot itself has been printed offline. 


L. ROOTS PROGRAM 
The ROOTS program is used to determine the roots of 
the characteristic equation for a general ARIMA model. The 


mOGeS are useful for testing for stationarity and for 
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determining the Pann of the forecast function. ROOTS uses 
the IMSL subroutine ZPOLR to calculate the roots; it requires 
as input the number of AR parameters in undifferenced form, - 
and the actual values of these parameters, as determined 

from the ESTIMATE or WMARQRDT programs. The roots are 
printed out at the terminal. If desired, roots can also 

be obtained for the MA polynomial, simply by substituting 


the appropriate MA values for the AR values as input. 


M. GENERATE AND SIMULATE PROGRAMS 

Because of their similarities, the GENERATE and SIMULATE 
programs are described together. The GENERATE program per- 
mits the user to generate a time series from any non-seasonal 
ARIMA model he specifies. The user must identify the model 
and give values for the model parameters and starting con- 
ditions; a random number seed must also be input. The 
program takes the specified model, generates random noise 
terms, and calculates as many values of the time series as 
desired. The GENERATE program can be useful for purposes 
of classroom instruction in the generation of a wide variety 


of time series examples for model identification exercises. 


tcuppose the model was ARIMA(1,1,0). The differenced form 
would be (1-9,B) Viz, = a + a,: Then, the undifferenced form 
would be written (1-(1+9,)B + 6B) 2, = of am ay , where the 
left-hand side was found by multiplying (1-9, 8B) by (i-B) = le 
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It could also be useful in the diagnostic phase of model 
checkout; a time series could be generated from the esti- 
mated model, and its properties compared with those of 

the original series. If large discrepancies occur in this 
comparison, this may be evidence that the model is inade- 
quate. The model output, the generated time series, is 
written in FILE FT02F001 as well as being printed offline. 
The GENERATE program uses IMSL subroutine FTGEN1. 

The SIMULATE program provides the capability of gen- 
erating any number of simulated time series. This program 
is useful for prediction of what might happen in the future; 
it also demonstrates that within a given model, the time 
series actually observed can vary considerably. This program 
uses GENERATE, but also requires as input the actual series 


1 and the number of simulated series the 


Sma FILE FTO2FOOL 
user wishes to generate. The series should be already 
transformed and differenced. Additionally, SIMULATE allows 
the user to select values of the original time series as 

| Starting values for the simulated series. The output con- 


sists of the simulated series printed at the terminal. The 


SIMULATE program uses IMSL subroutine FTGENI. 


mM. HELP OPTION 
The HELP option is a modification to the Time Series 


Editor CP/CMS program itself that allows the user to obtain 


Se desired, the user may choose to input his own starting 
values, and not use the actual time series values. 
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information about Editor programs after the initial 


introductory phase has been completed. 


O. SUMMARY 

This chapter has provided brief descriptions of the 
options available for analysis of time series data in the 
Time Series Editor. The following two chapters contain 
detailed examples of time series analysis for both a non- 
seasonal and a seasonal series; the accompanying terminal 
sessions and computer output for analysis of these series 
are contained in full in Appendices B and C. Appendix A 
to this report contains a User's Guide for the Tie Series 
Editor, while Appendix D contains a complete listing of 
the Time Series Editor and resident programs. Appendix E 
contains listings of the data sets used in the Hei sas o£ 


the next two chapters. 
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V. EXAMPLE NON-SEASONAL TIME SERIES ANALYSIS 


This chapter provides a description of a complete analy- 
sis of a non-seasonal time series using the Time Series 
Editor. The time series analyzed is Series C (Chemical 
Process Temperature Readings), taken from Box and Jenkins 
[Ref. 4, p. 528]. This time series was selected since it 
is fully analyzed and discussed in Reference 4, so that the 
interested user can compare for himself the results pre- 
sented there with the results given by the Time Series 
Editor. After reading this overview of the analysis, the 
reader can follow the actual user session and associated 
user session Editor program output contained in Appendix B 
Eo this report. 

mae first step in the analysis of series C is to plot 
the series uSing program PLOT. The plot reveals rather wide 
fluctuations in the series, but not the sort of explosive 
nonstationary behavior that would render an attempt to model 
eeureless. The plot also reveals that the time series has 
a large amount of momentum (movements of the series tend to 
resist changes of direction over time). This is character- 
istic of ARIMA models with two or more differences. 

As a second step, the autocorrelations (autos) and 
Partial autocorrelations (pautos), mean, and variance of 
the series were estimated using AUTO. The plot shows that 


the autos dampen out slowly in a nearly linear fashion. 


84 


Gl6yY [AMA Bale 


dren 4s Io nelsqe7O8eS a BEDE q 


ze tee ealt Lance 


tvisnha eeirtee 


| 


} 


This is an indication that the process is nonstationary 

and that one or more differences are required to make it 
stationary. The pauto plot is not informative when the 

autos fail to dampen out rapidly. 

As suggested by the plots of the original series and 
its auto and pautos, the program DIFF was used to examine 
the series with both one and two non-seasonal differences; 
i.e., the series zy = Viz, and Zz = Vz, were derived. 

The plots suggest that two differences might be required to 
achieve stationarity in the series. However, one might be 
able to get by with a single difference and an AR parameter 
near unity, Since the autos of zy dampen out slowly. Thus, 


two possible model candidates are suggested: 


im ARIMA(1,1,0) : (obs $B) V"z, = eS SF avy 


and 
MaMa (6,2 2ivaew 2. = 6. + (1 - 6.8 - 9B )a, 


For purposes of estimation, the second model was extended 

to two moving average parameters. Such overfitting is often 

done to see if the estimated. moving average parameters turn 

Out to be near zero, thus confirming the tentative identification. 
The next step in the analysis process is to calculate 


the maximum likelihood estimates of the model parameters. 
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The program ESTIMATE was used to do this. The estimated 


parameters for the ARIMA(1,1,0) model were: 


$, = 0.8073 
89 = -0.006789 
ae =. 1090177633 


The autos and pautos of the residuals were calculated 
to test the model. The correlations should appear to be 
estimates of a pure white noise process if the model is 
adequate; the (1,1,0) model seemed to pass that test. The 
chi-square goodness-of-fit test results ina es = 28.88 with 
24 df and a signficance level of 0.2247. Thus, there is 
no strong evidence tc suggest that the (1,1,0) model is 
inadequate. The model can be written, in undifferenced 


FOrM : 


Gie=—1.3073B + 0.8073B°)z, ont 


The 8, is omitted here since it is nearly zero in value. 


The estimated parameters for the (0,2,2) model were 


o5 = 051378 

oD) = 40512296 

8 = -0.0026787 
ce) 

of =e 00189559 


As before, the correlation plots of the residuals fail to 
Suggest any inadequacy of the model. However, the chi-square 
goodness-of-fit test for this model yielded a ei = 36.75 


with 23 df and a significance level of 0.0346. Therefore, 
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if the ARIMA(0,2,2) model were a "correct" one, a chi-square 
value as large as 36.75 would occur by chance with a proba- 
batity of only 0.0346. Due to its simplicity (parsimony 

is a desirable feature in time series models) and better 
eit, the ARIMA (1,1,0) model was: chosen as the "best" of 

the two alternatives originally selected. 

The ARIMA(1,1,0) model calculated was then used to make 
forecasts and confidence limits for these forecasts, using 
program FORECAST. The forecasts were made for 25 periods 
into the future, with plot origin 200 and forecast origin 
220. The length of the series C itself is 226; with one 
difference, it becomes 225. The plot in Appendix B shows 
the forecasted values, along with the confidence limits at 
level of significance 0.10; the plot shows that the width 
of the confidence limits increases substantially as the lag 
gets large. 

Finally the program ROOTS was used to calculate the 
roots of the characteristic equation for the ARIMA(1,1,0) 
model, which can be written: 


T= 1.80733: 0.8073B- = 0 


ine roots calculated by ROOTS were 1.239 and 1.0. Also, 

the program SIMULATE was used to generate a simulated series 
from the ARIMA(1,1,0) model, using the last two values of 
the actual series as starting values. The entire session 
lasted about two hours, including checkout of plots and 


consumed less than a minute of CPU time. 
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VI. EXAMPLE SEASONAL TIME SERIES ANALYSIS 


This chapter provides a description of a complete analy- 
sis of a seasonal time series using the Box-Jenkins tech- 
nique through the Time Series Editor. The time series 
analyzed is Series G (International Airline Passengers: 
Monthly Totals (Thousands of Passengers) January 1949 - 
December 1960), taken from Box and Jenkins [Ref. 4, p. 531]. 
This time series was selected since it was fully analyzed 
and discussed in Chapter 9 of Reference 4; again, the 
interested reader can compare for himself the results pre- 
sented by Box and Jenkins with those given by the Time 
Series Editor programs. After reading this overview of the 


analysis, the reader can follow the actual user session and 


associated Editor program output contained in Appendix C 


Borthis report. 

The plot shows both annual seasonality and an increasing 
trend. This suggests that a seasonal model with differencing 
may be needed. The autocorrelations (autos) and partial 
autocorrelations (pautos) were then calculated and plotted 
using AUTO. The plot of the autos shows very slow damping 
Out with peaks at intervals of 12, indicating strong non- 
stationarity and a seasonal period of 12 months. Again, 
here the pauto plot is not informative, due to the slow 
damping out of the autos. As it is sometimes useful in 


dealing with seasonal models, a natural log transform was 
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made on the series G data using program TRANS. The plot 
of the logged data shows little change in non-stationarity, 
and the autos and pautos of the logged data still indicate 
seasonality and non-stationarity. 

Next, differencing to achieve stationarity was tried, 
using program DIFF; both one seasonal and one non-seasonal 


difference of the logged data was taken. The autos and 


if 
12 


rapidly damping out. They also exhibit sharp peaks at the 


pautos of the data viy (ln z,) appear much improved, 
periods of seasonality, indicating the need for seasonal 
parameterization in the model. Since positive identification 
of the model type from the auto and pauto plots was not 


possible, two were postulated as close candidates: 


1. ARIMA ie) = ese , written: 


ales: Ke ey mn 12 
V Vy> (in z,) = 8G + (1 8,8) (1 ©, 128 a, ; 
and 
2 : 
oo ARTMAS (15.1;1) x (0,1;1) > Wratten: 
= 151 Ys 2 Es 12 
ell 8,B)V Vy5 (in Zz.) = 86 ei CE: 9,8) (1 Tae )a,- 


The next step in the analysis was to calculate the 


Parameter estimates for both models. YESTSEAS was used to 
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obtain initial parameter estimates. For the model 


ARIMA Rope) = (0,41) ~~, YESTSEAS yielded starting values: 


JE ey = 0.390425 


“a 


01 12(s) 02534397 . 


For this model, calculation of final parameter estimates 
in WMARQRDT using the starting values above yielded the 


following in seven iterations: 


ey =5,0437 7152 
ay waiwe- 282287 
6 = 0.000291 ; 
Oo 


the chi-square statistic for residual lack of fit of the 
model was ve =4292713571,.waithadiv=)38.andsa.probability 
of exceeding the eG value of .829. 

For the ARIMA (1,1,1) x (0,1,1)1* model YESTSEAS calcu- 


lated starting values: 


0.112699 


fap) 
— 
Qn 
~~ 
HT} 


0.490883 
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for this model, WMARQRDT yielded the following parameter 


values in ten iterations: 


$, = 0.146667 
a2 = 0.509816 
yee CRREEEE 


The chi-square statistic for residual lack of fit of the 
model was ae = 28.302543, with df = 37 and a probability 
of exceeding the Si value of .847. 

Next, the plots of autos and pautos of the residual for 
both models were examined; both models have autos and pautos 
that dampen out rapidly, and essentially appear to be 
"white noise", with the exception of a small peak at auto 
and pauto of lag 12, corresponding to the seasonal period. 
The chi-square statistics indicate that both models have 
"about" the same goodness of fit to the data. Consequently, 
both models were used to make forecasts of the series. 

Program FORECAST was used to forecast series G into the 
eaeure £or 30 time periods, with a forecast origin of 131; 
the plot origin selected was index 100. As expected, both 


models do about the same job of forecasting the series. 


Oa 


nn ‘oe 


The significance level for the confidence limits for both 
plots was 90%. 

As a final check of both models, the program XSUMSQ 
was used to calculate the sum of squared residuals for 


models. The results were as follows: 


67) = 0.281929 


all craton Js UilemlS Ws We mala 


De For ARIMA (a as de) x (G) elas tA 


A “ “a 


iS) (55 Oe ee = 0.181643 . 


It can be seen here that both models are nearly equally 
"good", and either could be used with a reasonable amount 


of confidence for forecasting series G. 
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VII. SUMMARY 


There has been an increasing need in military applica-. 
tions, as well as in the general business area, for a 
computer-aided time series analysis and forecasting capability 
that is both powerful and easy to use. Current military- 
associated functions in which time series techniques are 
valuable include material inventory Management, vecruiting, 
personnel management, budget analysis, force level projections, 
and perhaps even short-term weather prediction. 

The Time Series Editor that has been described in this 
report provides a unified collection of programs that greatly 
facilitate the conduct of a complete time series analysis 
and forecasting evolution using Box-Jenkins methodology. 

The routines guide the user from the data input stage, 
through the data analysis and parameter estimation stages, 
and finally through the model building, diagnostic checking 
and forecasting stages. With its simple input requirements, 
the Editor can be quickly mastered by even a beginning 
computer user, who has a basic understanding of Box-Jenkins 
methodology. 

The Box-Jenkins methodology has been described broadly 
in Chapter II and descriptions of the programs resident 
in the Time Series Editor have been given in Chapter IV. 
Example analyses of both non-seasonal and seasonal time 


series are provided in Chapters V and VI, with complete 


S)) 


user sessions and output packages for these analyses 
included as Appendices B and C. The Appendices also con- 
tain a Time Series Editor User's Guide, as well as complete 
program listings. 

Although the Time Series Editor in its present form 
provides an excellent capability for interactive time 
series analysis using the Box-Jenkins methodology, there 
are additions to the program set that might be made in the 
future which would improve the Editor's capability and 
utility. It is first recommended that the model diagnostic 
testing capability be extended to include a periodogram 
analysis, and/or other tests related to the spectral analy- 
sis of time series [see Ref. 4, p. 294]. Another useful 
addition might be to include in the Time Series Editor a 
package capable of improving the ease of entering data 


and simplifying the input error correction problem. 


94 


senpvirans Sacl}y se 


noe ' & ee>to0sqqA ser 


bom 


- 
, =>! tind 


APPENDIX A 


USER'S GUIDE TO THE TIME SERIES EDITOR 


The Time Series Editor is a collection of FORTRAN 
programs driven by a control program called TIMESER EXEC 
in the CP/CMS executive language that has been designed 
specifically for the analysis and forecasting of time 
series data using the general Box-Jenkins methodology. 
This Guide contains information essential for the user to 
access the Time Series Editor, enter time series data, 


build a model and obtain the desired output. 


I. DATA INPUT 

The Time Series Editor requires that time series data 
be entered into the sequential FORTRAN input/output file 
named FILE FTO2F001. This can be done either online or 
via cards read offline. The following diagram illustrates 
the proper card deck arrangement. When your deck is ready 
for input, give it to a computer center system operator 
to direct it to your disc space. 

If the user has a data deck already punched up ina 
format other than 5F15.6, he may enter the series as above 
into FILE FTO3F001 (without the length of series card) and 
use the Editor program ZFORMAT to transform it into the 
FILE FTO2F001 in the proper format, without destroying his 


Original file. This program is described in the next section. 
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II. TABLE OF OPTIONS 
This Table provides the user with the basic information 
necessary to understand the data requirements, functions 


and output options for each program in the Editor. 
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Tit. THE BASIC USER SESSION 

To use the Time Series Editor, the user must log into 
CMS, get into CP, link to the disc storage area where the 
Time Series Editor resides, reimplement CMS, log into the 
general user and Time Series Editor disc areas, and enter 
the TIMESER routine. This section will provide explicit 
guidelines to enable the user to perform the above steps 
on the NPS CP/CMS system. Commands marked with an asterisk 
(*) are those actually entered on the terminal by the user 
(the asterisk itself is omitted). Those without an aster- 
isk and those written in all capital letters are system 
responses at the terminal. Numbered sentences are comments, 
which will not appear during an actual user session. The 
instructions and system responses assure the user is on an 
IBM 2741 Input/Output Terminal. Some minor modifications 


may be necessary if other terminals are used. 


1. Turn the terminal on, depress the BREAK key, and wait 
for the system to respond: 


CP-67 online xd.65 gsyosu 


2. Depress the ATTN key. The roll bar will advance and 
the keyboard will unlock. MThen enter: 


*login aaaapbb 450k 


3. aaaa is the user's identification number, and nn is the 
terminal number (usually written on the terminal). 
For example, if the user's ID number is 1621 and the 
terminal number is 44, the input would be: 
Hogin, 1621p44 450k. The addition of 450 k to the normal 
login command is necessary to execute the program 
WMARQRDT in the Editor; for users not planning to 
execute this program during a session, this addition 
is not necessary. 
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The system will respond with the statement: 
ENTER PASSWORD: 


The user then enters his password, or the general users 
password npg; 


*password 
The system will then respond: 


ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER 
COST CENTER CODE: 


The user then enters: 
*gggghhhh 


gggg is the assigned project number, and hhhh is the 
user's section designator or the faculty code. 


The system will respond with the message of the day, 
such as: 


P7CMS HOURS ... 0930=2200 (MON-THURS) ... 0930-1800 (FRI) 


“OUTPUT RETAINED 5 DAYS 


Cms Version 3.25 

At this point, the user is in CMS. He must then get into 
CP; this can be done by hitting the ATTN key. The system 
will then respond: 

CP 


The user must then link to the TIME SERIES EDITOR; 
this is accomplished by entering: 


Sienk 2069p 191 192 

The system will respond with: 

ENTER PASSWORD: 

The password (read only) to enter the Editor is: 
*timser 

The system then responds: 


SET TO READ ONLY 
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The user now implements CMS by: 
*ipl cms | 

The system will respond: 

CMS Version 3.25 


Now the user must log into both the general user and 
the Time Series Editor area by entering: 


*login 191 

The system ae respond with a message such as: 
R; 

The user then enters the command: 

‘etogin 192 t,p 


The system will respond: 


mt l92). R/O 


’ 


The user can then enter the Time Series Editor 
(guided version) by entering the command: 


*timeser 
The system will respond: 
mach 2 SECONDS EXECUTION TIME IS INDICATED BY * 


YOU HAVE ENTERED THE TIME SERIES EDITOR 


PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 


ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 
ENTER NUMERICAL VALUES VIA FORTRAN FORMAT., 


TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THROUGH N. TYPE FLOATING VALUES WITH DECIMAL FOR 
ALL OTHERS. 

be YOU WANT A. LIST OF THE OPTIONS? 

The user is then on his own, guided by the Exec routine. 
See the notes that appear at the end of this guide for 
additional information. Eventually the user will be 
asked: 


DO YOU WANT TO TRY AGAIN? 


LOW 


Dat | 


i 


__— 


24. 


25. 


26. 


i 


If a yes response is given, another sequence will begin; 
if the response is no, the user will be taken out of 
the Time Series Editor environment and returned to CMS. 
The system response will be: 


CONTROL RETURNED TO CMS 
R; 


The user can then log out of CMS by typing: 
*cp logout 
The system will respond with: 


CONNECT= 00:08:02 VIRTCPU= 000:07.98 TOTCPU= 000.10.94 
LOGOUT AT 14.22.04 on 10/16/78 


The user should then turn off his terminal and tear off 
the output from his session. 


The more experienced user can dispense with the 


"welcome aboard" section of the Time Series Editor and get 


right down to business by using the shortened version of the 


Editor. This shortened version may be entered by linking in 


the normal way, and then entering the Editor by typing the 


COMMAND 


*timeser s (asterisk omitted). 


The system will immediately respond: 


ENTER LETTER FOR OPTION YOU WANT. 


The session inside the Editor then begins. 
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IV. BRIEF SAMPLE USER SESSION 


A brief sample user session is given below; it includes 


copies of the offline output generated during the session. 
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repeat login nur@pné@ 


login 1621p44 450k 
ENTER PASSWORD: 


ENTER 4-DIGIT, PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0444r172 ° 

READY AT 17.24.38 ON 09/16/78 

CMS Version 3.25 


stat 

P (191): 29 FILES; 241 REC IN USE, 55 LEFT (of 296), 81% FULL (2 CYL) 
R; : 

cpqf 

FILES:- NO RDR, NO PRT, NO PUN 

R; 

CP 

link 2069p 191 192 

ENTER PASSWORD: 


SET TO READ ONLY 


ipl cms 
CMS Version 3.25 


login 191 

R; 

login 192 t,p 

T (192) R/O 

R; 

timeser 

EACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 


YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTeR ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 


TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THRU N. TYPE FLOATING VALUES WITH DECIMAL FOR ALL OTHERS. 


DO YOU WANT A LIST OF THE OPTIONS? 


y 

OPTION DESCRIPTION 

GENERATE ----- GENERATE ANY ARIMA TIME SERIES 

[\ ONG) CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 
BEOLe=——-=———=—=PLOT A TIME SERIES 

ESTIMATE ----- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 
OME) OS DIFFERENCE A TIME SERIES 

FORECAST ----- FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 
TRANS ------- “TRANSFORMS VALUES OF A TIME SERIES 

ROOTS -------- DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 
ZFORMAT ------ ALTER DATA FILE TO FORMAT 5F15.6 

CMSWORK ------ PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE ----- SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS ----- CALCULATE INITIAL SEASONAL PARAMETERS 

wHARQRDT ----- MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ ------- CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 


WOULD YOU LIKE MORE INFO? 
yi 
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ENTER OPTION YOU WANT INFO ABOUT. 


AUTO -- 


WHICH MUST RESIDE IN FILE FTO2FO01. 
FTAUTO IN THE IMSL LIBRARY. 


DSRS TS THIS PROGRAM CALCULATES AUTOCORRELATIONS, PARTIAL 
AUTOCORRELATIONS, THE MEAN AND THE VARIANCE FOR A GIVEN TIME SERIES 
THE PROGRAM USES 


PLOTTED OFFLINE. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


n 


DO YOU WANT TO TRY A SESSION? 


¥ 


ENTER LETTER FOR OPTION YOU WANT. 


Cc 


ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
| WHEN FINISHED TYPE: &GOTO -QUES 
listf * f£t02f001 


FILENAME FILETYPE MODE NO.REC. DATE 


SERG 
SERC 


LNSERG 


FILE 


FTO2FOOl Pl 3 9/16 
FTO2FO0O1l Pl 5 9/16 
FTO2FOO1 Pl 3 9/16 
FTO2FO01 Pl 1 9/16 


erase file f£t02f001 
alter serc £t02f001 pl file ft02f001 pl 


stat 


P (191): 28 FILES; 240 REC IN USE, 
&goto -ques 
DO YOU WANT TO GO AGAIN? 


¥ 


ENTER LETTER FOR OPTION YOU WANT. 


a 


IS YOUR DATA IN FILE FTO2F001? 


Mes. 
EXECUTION BEGINS... 


AUTOCORRELATIONS 

Ges7a 0.944 0.902 0.854 0.802 
Deaee 0.413 0.359 0.305 0.253 
=0.052 -0.101 -0.151 -0.200 -0.248 

PARTIAL AUTOCORRELATIONS 
0.978 -0.260 -0.157 -0.093 -0.058 
~=0.036 -0.041 -0.038 -0.024 -0.037 
-0.034 -0.061 -0.079 -0.048 -0.037 
MEAN= 22.9739 VARIANCE 


ENTER TITLE FOR PLOTS. 

autos and pautos for series c data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 
DO YOU WANT TO GO AGAIN? 


¥ 


ENTER LETTER FOR OPTION YOU WANT. 


56 LEFT (of 296), 


0.748 0.692 
0.201 0.150 


0.635 
0.098 


THE AUTOCORRELATIONS AND PAUTOS CAN BE 


81% FULL (2 -CYL) 


0.579 0.323 
0.047 -0.003 


-0.045 -0.012 -0.038 -0.022 -0.010 
-0.027 -0.032 -0.070 -0.048 -0.024 


4.22273 


c 

. ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

’ offline print file f£t02£001 

alter file £t02£001 pl serc ft02f001 pl 

&Goto -ques 

DO YOU WANT TO GO AGAIN? 


n 
CONTROL RETURNED TO CMS 


R; 
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V. PROBLEM CONTROL NOTES 

This section will cover corrective measures that can be 
taken when things fail to go as expected while in the 
TIMESER environment. 

a. A typing error in the CMS environment can be 
corrected by typing the @ character as many times as is 
required to back up and then type the correct values. For 
example, if the user typed timesre, the user could correct 
the mistake by typing two @ signs, followed by the correct 
epetiing er, as follows: timesre@@re. An entire line can 
be deleted by typing the ¢ character (or [ on some terminals). 

b. When working with TIMESER executive programs, the 
user should exercise care before hitting the return key (or 
control s on some terminals). If an input value is required 
and the return key is hit before the proper response is 
entered, the user will be likely to get thrown out of the 
editor and have to begin again. In most cases, errors can 
be corrected only before the return key is struck (in some 
programs you get a second chance for input). 

c. Particular care should be taken for integer value 
input, which must be right justified in the format field. 
The editor will advise the user in all cases where the 
integer format is other than Il. 

d. If for any reason the user finds himself in a debug 
or error condition (caused by erroneous data, or a "blowup" 


in one of the non-linear optimization routines usually 


16 


oe 
. 


seven Die 


os tint 


caused by very poor initial input values), the following 
procedure will get the user back into the normal TIMESER 
environment: 
(1) depress the ATTN key twice; this gets 
the user into CP; hit the ATTN key again, 
and then type kx to kill the execution. 
(2)) re-ipl CMS, and login 191 and then login 
1O2 ep. a 
(3) then the user can type TIMESER or TIMESER S, 
and return to the TIMESER environment 
On the next page is a sample user session where an error 
cauSing a debug condition has occurred. 
Here the user executed a program that required data 
ea PrLE FTO2ZFO01lL, and the FILE did not exist. As the 
example shows, recovery is quick. Simply hit the break 


Baceon tO get into CP, login 191, login 192 t,p, and then 


type timeser s to return immediately to the Editor environment. 
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APPENDIX B 
KEYBOARD SESSION AND OUTPUT FROM ANALYSIS OF 
TIME SERIES C (NON-SEASONAL) 

This appendix contains the complete terminal keyboard 
session and all offline output from the analysis of Box 
and Jenkins' [Ref. 4] time series C (Chemical process 
temperature readings) using Cnet ine Series Editor, This 
appendix is intended to supplement the discussion of the 
analysis of series C discussed in Chapter V of this report. 

The terminal listings are presented first, followed by 


the computer offline output. 
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‘ti sucruo enna 


repeat login nur@pne@ 


login 1621p44 400k 
ENTER PASSWORD: 


ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0444r172 

READY AT 14.23.02 ON 09/16/78 

CMS Version 3.25 


listf * £t02f£001 
FILENAME FILETYPE MODE NO.REC. DATE 


FILE FTO2FOO1l Pl 3 9/15 
SERG FTO2FO01l Pl 3 9/16 
SERC FTO2F001 Pl 5 9/16 

R; 

erase file £t02f001 

Rj 

alter serc ft02f001 pl file £t02f£001 pl 
R; : 

CP 


link 2069p 191 192 
ENTER PASSWORD: 


SET TO READ ONLY — 


ipl cms 
CMS Version 3.25 


login 191 

R; 

login 192 t,p 
T (192) R/O 


4 


.er 
waCH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 


YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 


TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THRU N. TYPE FLOATING VALUES WITH DECIMAL FOR ALL OTHERS. 


DO YOU WANT A LIST OF THE OPTIONS? 


Y 

OPTION DESCRIPTION 

GENERATE ----- GENERATE ANY ARIMA TIME SERIES 

AUTO -<-----~--- CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 
PLOT --------- PLOT A TIME SERIES 

ESTIMATE ----- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 
DIFF --------- DIFFERENCE A TIME SERIES 

FORECAST ----- FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 
TRANS 9-------- TRANSFORMS VALUES OF A TIME SERIES 

ROOTS -------- DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 
ZFORMAT ------ ALTER DATA FILE TO FORMAT 5F15.6 

CMSWORK =----- PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE ----- SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS ----- CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT ----- MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ ------- CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 


WOULD YOU LIKE MORE INFO? 
Y 
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ENTER OPTION YOU WANT INFO ABOUT. 
a 


AUTO --------THIS PROGRAM CALCULATES AUTOCORRELATIONS, PARTIAL 
AUTOCORRELATIONS, THE MEAN AND THE VARIANCE FOR A GIVEN TIME SERIES 
WHICH MUST RESIDE IN FILE FTO2F001. THE PROGRAM USES 

FTAUTO IN THE IMSL LIBRARY. THE AUTOCORRELATIONS AND PAUTOS CAN BE 
PLOTTED OFFLINE. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


Y 

ENTER OPTION YOU WANT INFO ABOUT. 

d 

POSS! Oe THIS PROGRAM TAKES A SERIES IN FILE FTO2F001, ASKS IF 


IT IS SEASONAL, AND THEN PERFORMS NONSEASONAL AND/OR SEASONAL 
DIFFERENCING, AS REQUESTED. THE ORIGINAL DATA REMAINS IN FILE 
FTO2F001, WHILE THE DIFFERENCED SERIES IS PUT IN FILE FTO3FO001l. 
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DO YOU WANT INFO ABOUT ANOTHER OPTION? 
n 
DO YOU WANT TO TRY A SESSION? 


Y 
ENTER LETTER FOR OPTION YOU WANT. 


a 
‘IS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


On978 0.944 0.902 0.854 0.802 0.748 0.692 0.635 0.579 0.523 
Geaes 02413 0.359 0.305 0.253 0.201 0.150 0.098 0.047 ~0.003 
=0.052 —-0.101 -0.151 -0.200 -0.248 


PARTIAL AUTOCORRELATIONS 


0.978 -0.260 0.157 -0.093 -0.058 -0.045 -0.012 -0.038 -0.022 -0.010 
0.036 -0.041 -0.038 -0.024 -0.037 -0.027 -0.032 -0.070 -0.048 -0.024 
-0.034 -0.061 -0.079 -0.048 -0.037 


MEAN= 22.9739 VARIANCE = 4,22273 


ENTER TITLE FOR PLOTS. 

autos and pautos of series c data / undifferenced 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 

*DO YOU WANT TO GO AGAIN? : 

YY, 
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ENTER LETTER FOR OPTION YOU WANT. 
d ; ‘ 
IS YOUR DATA IN FILE FTO2F001? 


Y, 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

1 

*DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


y 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.805 0.653 0.526 0.442 0.380 0.318 0.262 0.186 0.139 
-0.097 0.094 0.074 0.073 0.070 0.072 0.089 0.048 0.041 
G.044 0,048 0.000 -0.052 -0.097 


PARTIAL AUTOCORRELATIONS 
0.805 0.010 -0.007 0.051 0.027 -0.019 -0.013 -0.080 0.020 
-0.137 0.094 -0.027 0.034 0.006 0.012 0.043 -0.121 0.060 
0.038 -0.022 -0.116 -0.076 -0.008 
MEAN=-.346667E-01 VARIANCE = 0.531982E-01 
ENTER TITLE FOR PLOTS. 
autos and pautos of series c data / one nonseasonal difference 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 


Y 

ENTER LETTER FOR OPTION YOU WANT. 

d 

IS YOUR DATA IN FILE FTO2F001? ~ 


yy, 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 
n 
ENTER NUMBER OF NONSEASONAL DIFFERENCES. 
2 


DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


MY 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


-0.079 -0.065 -0.122 -0.063 0.013 -0.018 0.049 -0.052 -0.124 
soei22 0.072 -0.077 0.029 -0.011 -0.058 0.171 -0.101 -0.013 
=0.007 0.135 0.014 -0.013 -0.135 


PARTIAL AUTOCORRELATIONS 
-0.079 -0.072 -0.135 -0.094 -0.022 -0.051 0.022 -0.060 -0.145 
-0.143 0.023 -0.089 -0.005 -0.039 -0.074 0.130 -0.105 -0.014 
0.004 0.081 0.056 -0.041 -0.099 
MEAN=~.267866E-02 VARIANCE = 0.198143E-01 
ENTER TITLE FOR PLOTS. 
autos and pautos of series c data / two nonseasonal differences 


YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 
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DO YOU WANT TO GO AGAIN? 
y , 


* 


ENTER LETTER FOR OPTION YOU WANT. 


c 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE:, &GOTO -QUES 

alter file f£t02£001 pl serc ft02f001 pl 
alter file £t03£001 pl file ft02£001 pl 
&goto -ques 


DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
e 

IS YOUR DATA IN FILE FTO2F001? 

y 

EXECUTION BEGINS... 


0 
ENTER NUMBER OF MA PARAMETERS. 
2 


LENGTH OF TIME SERIES = 224 
0 AR PARAMETERS 2 MA PARAMETERS 


MA PARAMETERS ARE: 
THETA(1)= 0.1378 


THETA(2)= 0.1296 


ENTER NUMBER OF AR PARAMETERS (EXCLUDE DIFFERENCES). 


MA CONSTANT= -0.26787E-02 WHITE NOISE VARIANCE= 0.189559E-01 


CHI SQUARE RESIDUAL LACK OF FIT VALUE = 36.75 


SIGNIFICANCE = 0.0346 


DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 


Y 
*EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.019 0.034 -0.134 -0.093 -0.012 -0.043 0.005 
=O-a42 0.069 -0.098 0.025 0.004 -0.055 0.157 


0.009 0.133 0.005 -0.015 -0.133 


PARTIAL AUTOCORRELATIONS 


02019 0.033 -0.135 -0.090 0.001 <-0.055 -0.019 
=-0.166 0.017 -0.104 -0.013 -0.025 -0.089 0.123 


0.024 0.069 0.024 -0.069 -0.102 


MEAN=0.236684E-02 VARIANCE = 0.189503E-01 


ENTER TITLE FOR PLOTS. 


DF = 23 


-0.:063 -0.125 0.088 
“0.090 0.005 -0.006 


-0.071 -0.143 0.087 
-0.123 -0.021 -0.025 


autos and pautos of residuals / arima(0,2,2) model of series c data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 


PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
¥ 


£23 


jo'taa see 
" 4 FI TMD 72 
. ’ 1 oF 


a 
vitae On 
7 , 


4’ (0.9 7 
0- eat. a 
0 0 EE E.@ 


A Ac Tae 
-c1.0~\08 Rae 
bo! .G=— Tia 
60.0 26 


= — * 


Tal? 80% 
cvubiest 2o oe i. 
it OTA Gea 

U Ob17 SOG Ge 


° 
"wIeoA O38 OF 


ENTER LETTER FOR OPTION YOU WANT. 

c 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

erase file £t02f£001 

stat 

P (191): 27 FILES; 265 REC IN USE, 31 LEFT (of 296), 90% FULL (2 CYL) 
erase file ft08f001 

alter serc f£t02£001 pl file ft02f£001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

Y 


ENTER LETTER FOR OPTION YOU WANT. 
d 
*IS YOUR DATA IN FILE FTO2F001? 


bf 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 

bs 3 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 


1 
DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


n 
DO YOU WANT TO GO AGAIN? 
y 


ENTER LETTER FOR OPTION YOU WANT. 


c 

*ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

alter file £t02£001 pl serc ft02f001 pl 

alter file £t03f001 pl file ft02f001 pl 

&goto -ques 

DO YOU WANT TO GO AGAIN? 

sf 


ENTER LETTER FOR OPTION YOU WANT. 
e 


IS YOUR DATA IN FILE FTO2F001? 


yi 
EXECUTION BEGINS... 
ENTER NUMBER OF AR PARAMETERS (EXCLUDE DIFFERENCES). 
1 
ENTER NUMBER OF MA PARAMETERS. 
0 


* LENGTH OF TIME SERIES = 225 
1 AR PARAMETERS 0 MA PARAMETERS 


AR PARAMETERS ARE: 
PHI(1)= 0.8073 


MA CONSTANT= -0.66789E-02 WHITE NOISE VARIANCE= 0.177633E-01 


(CHI SQUARE RESIDUAL LACK OF FIT VALUE = 28. 88 DF = 24 
SIGNIFICANCE = 0.2247 - Yon 


- DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 


y 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.016 0.014 -0.049 -0.005 0.061 0.024 0.078 =-0.018 -0.087 0.130 


0.093 0.084 -0.055 0.043 0.007% -0.040 0.169 -0.079 =-0.001 =-0.010 
0.000 0.129 0.014 -0.013 -0.127 
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PARTIAL AUTOCORRELATIONS 


0.016 0.014 -0.050 -0.004 0.063 0.019 0.075 -0.015 -0.087 0.14° 
-0.104 0.070 -0.048 0.042 0.003 -0.030 0.158 -0.088 0.013 -0.025 
0.030 0.086 0.039 -0.062 -0.101 


MEAN=-, 252517E-02 VARIANCE = 0.177569E-01 


ENTER TITLE FOR PLOTS. 
autos and pautos of residuals / arima(1,1,0) model of series c data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
y, 


x 


ENTER LETTER FOR OPTION YOU WANT. 
Cc 


ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

erase file ft02f001 

alter serc ft02f001 pl file f£t02f001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 


ENTER LETTER FOR OPTION YOU WANT. 
6 

IS YOUR DATA IN FILE FTO2F001? 

yi; 

EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 
n 


ENTER NONSEASONAL DATA AS REQUESTED: 
ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

I! . 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 
1 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 
0 


NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, THETAO. 


-.00667 

ENTER NON-SEASONAL AR PARAMETER PHI(1). 

0.8073 

ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT I2, 
25 

ENTER INDEX FOR PLOT ORIGIN VIA FORMAT I3. 

200 

ENTER INDEX POR FORECAST ORIGIN VIA FORMAT I3. 
220 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 
0.10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
n : 


DO YOU WANT BASIC OUTPUT AT THE TERMINAL? 
Y 
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THE LAST 10 WVEC VALUES: 


0.222000E 02 0.218000E 


0.197000E 


02 


0.193000E 


THE 25 FORECAST VALUES: 
f 


0.203896E 
0.191821E 
0.186541E 
0.183591E 
0.181439E 


02 
02 
02 
02 
02 


THE 25 UPPER 


0.206094E 
0.206551E 
0.212721E 
0.219221E 
0.224989E 


ALPHA FOR THE 


02 
02 
02 
02 
02 


0.200S517E 
0.190434E 
0.185837E 
0.183121E 
0.181050E 


02 
02 


02 
02 
02 
02 
02 


0.213000E 
0.191000E 


0.197722E 
0.189248E 
0.185203E 
0.182676E 
0.180669E 


02 
02 


02 
02 
02 
02 
02 


FORECAST CONFIDENCE LIMITS: 


0.205056E 
0.207621E 
0.214056E 
0.220442E 
0.226043E 


02 
02 
02 
02 
02 


0.204779E 
0.208812E 
0.215381E 
0,221630E 
0.227064E 


02 
02 
02 
02 
02 


CONFIDENCE LIMITS iS: 0.100 


DO YOU WANT TO GO AGAIN? 


Y 


ENTER LETTER FOR OPTION YOU WANT. 


r 
EXECUTION BEGINS... 
SOLVES FOR ROOTS OF CHARACTERISTIC EQUATION. 


ENTER NO. OF AR PARAMETERS (UNDIFF. FORM). 


2 


ENTER AUTOREGRESSIVE PARAMETERS: 


ENTER PHI(1). 


1.8073 


ENTER PHI(2). 


-0.8073 


THE 2 ROOTS OF CHARACTERISTIC EQ. 


1.239 
1.000 


0.0 
0.0 


DO YOU WANT TO GO AGAIN? 


Y 


e2s9 
1.000 


ENTER LETTER FOR OPTION YOU WANT. 


d 


*EXECUTION BEGINS... 


IS YOUR TIME SERIES SEASONAL? 


n 


ENTER NUMBER OF NONSEASONAL DIFFERENCES. 
10} ° 


0.208000E 
0.190000E 


0.195399E 
0.188223E 
0.184623E 
0.182249E 
0.180295E 


0.205033E 
0.210079E 
0.216687E 
0.222783E 
0.228055E 


AND NORMS ARE: 


02 
02 


02 
02 
02 
02 
02 


02 
02 
02 
02 
02 


DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


n 


DO YOU WANT TO GO AGAIN? 


¥ 


ENTER LETTER FOR OPTION YOU WANT. 


Ss 


*EXECUTION BEGINS... 
lade NUMBER OF AR PARAMETERS (UNDIFFERENCED FORM) 


ENTER ESTIMATED AR PARAMETER PHI(1). 


1.8073 


ENTER ESTIMATED AR PARAMETER PHI(2). 


~0.8073 
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0.202000E 
0.188000E 


0.193456E 
0.187329E 
0.184089E 
0.181838E 
0.179926E 


0.205662E 
0.211390E 
0.217969E 
0.223903E 
0.229016E 


02 
02 


02 
02 
02 
02 
02 


02 
02 
O02". 
02 
02 
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AR PARAMETERS ARE 
ARE THESE OK? 


1.8073 


y ; 
ENTER NUMBER OF MA PARAMETERS 
0 


ENTER OVERALL MA CONSTANT -- 
-0.00667 


‘ENTER ESTIMATED WHITE NOISE VAR 
0.0177 


-0.8073 


DO YOU WANT STARTING VALUES TO BE THE LAST VALUES 


OF THE ACTUAL SERIES? 


Y 

START(1)= 19.000000 
START (2)= 18.799988 
STARTING’VALUE(S) OK? 


Y 
ENTER RANDOM NUMBER SEED (BETWEEN O AND 1). 


0.12345 


ENTER NUMBER OF VALUES YOU WANT TO SIMULATE VIA 13. 


50 


DO YOU WANT A PRINTOUT OF THE FIRST 20 VALUES? 


Wie. 
SIMULATED VALUES ARE: 


18.50224 18.33957 
17.85420 17.69398 
16.94655 16.68689 
16.31497 16.64629 


DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 
R; eas 


18.33694 
17.58267 
16.43227 
17.04015 
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18.16435 
1750722 
16:42499 
17.47137 


17.98804 
17.26422 
16.26595. 
17.62828 
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VARIANCE = 
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RESICLAL VALUES. 
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AUTCCCRRELATICNS 
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CRIGINAL ANC FORECASTED TIME SERIES. 
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SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS: 


FCRECAST CRIGIN:2220 


25 


MAXIMUM FCFECAST LEAD TIME: 


APPENDIX C 
KEYBOARD SESSION AND OUTPUT FROM THE ANALYSIS OF 
TIME SERIES G (SEASONAL) 

This appendix contains the complete terminal keyboard 
session and all offline output from the analysis of Box 
and Jenkins’ [Ref. 4] time series G (International airline 
passengers, monthly totals, Jan 1949 - pee 1960) using the 
Time Series Editor. This appendix is intended to supplement 
the discussion of the analysis of series G given in Chapter 
oeeot this report. 

The terminal listings are presented first, followed by 


the computer offline output. 
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32 SMrIT 


+ 


repeat login nuré@pn@ 


login 1621p44 400k 
ENTER PASSWORD: 


ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0444r172 

READY AT 15.12.34 ON 09/16/78 

CMS Version 3.25 


cP 
link 2069p 191 192 
ENTER PASSWORD: 


SET TO READ ONLY 


ipl cms 
CMS Version 3.25 


login 191 


R; 

login 192 t,p é 3 “ "OM ALAS 
T (192) R/O 

R; 


timeser 
EACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 


YOU HAVE ENTERED THE ous SERIES Bee x 


PLEASE RESPOND TO “BACH ‘QUERY WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LSTTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN ee. 

TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 

WITH I THRU N. TYPE FLOATING VALUES. WITH DECIMAL FOR ALL OTHERS. 


DO YOU WANT A LIST OF THE OPPIONS? . ‘ 

Y . 

OPTION DESCRIPTION 

GENERATE -~--- GENERATE ANY ARIMA TIMB SERIES | 

AUTO ----~---- -CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 
PLOT ~--------PLOT A TIME SERIES 

ESTIMATE ----- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 
DIFF --------- DIFFERENCE A TIME SERIES - 

FORECAST ----- FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 
TRANS ------>-TRANSFORMS VALUES OF A TIME SERIES 

ROOTS -------- DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 
ZFORMAT ------ ALTER DATA FILE TO FORMAT 5F15.% 

CMSWORK ------ PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE ----- SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS -----CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT ----- MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ ------- CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 


WOULD YOU LIKE MORE INFO?. | ° Sate 


y ‘ x oer a.’ 

ENTER OPTION YOU WANT INFO ABOUT. 

y 

YESTSEAS ----THIS PROGRAM CALCULATES INITIAL PARAMETER ESTIMATES 


FOR SEASONAL AND NONSEASONAL ARIMA MODELS, TO BE USED AS INPUTS TO 
THE WMARQRDT PROGRAM. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


y 
ENTER OPTION YOU WANT INFO ABOUT. 
Ww 


150. 


JJ ar 
a ,arveod Liv) 


2s 


£280 


WMARQRDT =---THIS PROGRAM CALCULATES NON-LINEAR LEAST SQUARES 
ESTIMATES OF BOX-JENKINS PARAMETERS FOR A GENERAL (SEASONAL 
OR NONSEASONAL) ARIMA MODEL. IT REQUIRES INITIAL PARAMETER 
ESTIMATES AS STARTING VALUES; THESE MAY BE CALCULATED USING 
PROGRAM YESTSEAS. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


y 
ENTER OPTION YOU WANT INFO ABOUT. 
x 


asUNsQ) =-—---- THIS PROGRAM CALCULATES THE SUM OF SQUARED RESIDUALS 
FOR ANY SET OF ARIMA (SEASONAL OR NONSEASONAL) PARAMETERS. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


n 
*DO YOU WANT TO TRY A SESSION? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 

Cc 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FPINISHED TYPE: &GOTO -QUES 

stat 

P (191): 26 FILES; 235 REC IN USE, 61 LEFT (of 296), 79% FULL (2 CYL) 
cp q £ 

FILES:- NO RDR, NO PRT, NO PUN 

listf * ft02f001 

FILENAME FILETYPE MODE NO.REC. DATE 

SERG FTO2FOO1l Pl 3 9/16 

SERC FTO2FO0O1 Pl eS) e710 

alter serg £t02f£001 pl file ft02f001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 


Y - 
ENTER LETTER FOR OPTION YOU WANT. 


P 
IS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 

ENTER TITLE FOR PLOT 

series g data / undifferenced 

TIME SERIES PLOTS HAVE BEEN PRINTED OFFLINE 
DO YOU WANT TO GO AGAIN? 
y 


ENTER LETTER FOR OPTION YOU WANT. 
a 
IS YOUR DATA IN FILE FTO2F001? 


¥ 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


@eeeeoearo 6€©60.807 0.753 0.714 0.682 0.663 0.656 0.671 0.703 
0.743 0.760 0.713 0.646 0.586 0.538 0.500 0.469 0.450 0.442 
0.457 0.482 0.517 0.532 0.494 

PARTIAL AUTOCORRELATIONS 
0.948 -0.229 0.038 0.094 0.074 0.008 0.126 0.090 0.232 0.166 


O.171 -0.135 -0.540 -0.027 0.091 0.025 0.032 0.077 0.048 -0.046 
0.046 -0.100 0.052 0.048 -0.163 
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pee 


/  SAMOTGE, A 


ro AoW 
ie 

‘4 7 at 

«es 

Q OTe 


ro6.0 
0 E228 
o 1F 8 


>€23.6 y 0.0 
feo ,o~ ‘ 7 7 
4>0.0 . 


MEAN= 280.299 VARIANCE = 14292.0 


ENTER TITLE FOR PLOTS. 

autos and pautos of series g data / undifferenced 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 


y 


ENTER LETTER FOR OPTION YOU WANT. 
t 
IS YOUR DATA IN FILE FTO2F001? 


y- 
EXECUTION BEGINS... 

DO YOU WANT-TO TRANSLATE THE ORIGIN: W=Z-SHIFT? 
n 

DO YOU WANT TO RESCALE THE VALUES: W=2*SCALE? 


n 

DO YOU WANT A LOG TRANSFORMATION? 

¥ 

TRANSFORMATION IS W(I)#LOG(SCALE*(Z(1I)-SHIFT)-FACTOR) WHERE: 


SCALE= 1.00000 SHIFT=0.0 FACTOR2=0.0 
YOUR ORIGINAL SERIES IS IN DATA FTO2F001 Pl. 
YOUR TRANSFORMED SERIES IS IN FILE FTO2F001 Pl. 


DO YOU WANT TO PLOT THE TRANSFORMED VALUES? 


iy 

EXECUTION BEGINS... 

* ENTER TITLE FOR PLOT 

plot of series g data / natural log transform 
TIME SERIES PLOTS HAVE BEEN PRINTED OFFLINE 


DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
a 
IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


Begoe 0.899 0.851 0.808 0.779 0.756 0.738 90.727 0.734 
aeaoe 0.762 0.717 0.663 0.618 0.576 0.544 0.519 90.501 
Oeage 0.506 0.517 0.520 0.484 - 


PARTIAL AUTOCORRELATIONS 


0.954 -0.118 0.054 0.024 0.116 0.044 0.038 0.100 0.204 
0.106 -0.042 -0.485 -0.034 0.042 -0.044 0.028 0.037 0.042 
0.073 -0.033 0.061 0.031 -0.194 a 


MEAN= 5.54218 VARIANCE = 0.193531 


ENTER TITLE FOR PLOTS. 

autos and pautos of series g data / natural log transform 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

Y 


$52 


0.744 
0.490 


0.064 
0.014 


org Sate aE 


‘TALUWART 
anew OF 
iT Tw) @ 
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ENTER LETTER FOR OPTION YOU WANT. 
d 
IS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


Y ENTER ORDER OF SEASONAL DIFFERENCING. 

= LENGTH OF SEASONAL PERIOD VIA 12. 

ss NUMBER OF NONSEASONAL DIFFERENCES. 

4 YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


y, 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


“0.341 0.105 -0.202 0.021 0.056 0.031 -0.056 -0.001 0.176 
0.064 -0.387 0.152 -0.058 0.150 -0.139 0.071 0.016 -0.011 
0.039 -0.091 0.223 -0.018 -0.100 


PARTIAL AUTOCORRELATIONS 


-0.341 -0.013 -0.193 -0.125 0.033 0.035 -0.060 ~0.020 0.226 
0.047 -0.339 -0.109 -0.077 -0.022 -0.140 0.026 0.115 -0.013 
0.132 -0.072 0.143 -0.067 -0.103 


MEAN=0.290920E-03 VARIANCE = 0.208604E-02 


ENTER TITLE FOR PLOTS. 

autos and pautos for series g s@data / 1 ns diff, 1 seas diff, 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 

c 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

alter file ft03f£001 pl save ft03f001 pl 
alter data f£t02f£001 pl serg ft02f001 pl 

stat 


-0.076 
-0.117 


0.043 
-0.167 


ln xform 


Bet29l): 29 FILES; 267 REC IN USE, 29 LEFT (of 296), 90% FULL (2 CYL) 


erase file £t08f£001 

ep qf 

FILES: = NO RDR, NO PRT, NO PUN 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 


ENTER LETTER FOR OPTION YOU WANT. 
Y 
IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 
oo LENGTH OF SEASON VIA FORMAT I2. 


— NUMBER OF NON-SEASONAL DIFFERENCES. 
hala NUMBER OF SEASONAL DIFFERENCES. 


_ NUMBER OF NON-SEASONAL AR PARAMETERS. 
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oosee : 
:\088093 = 


Ou wa 


on / oF) “ 
i ; 


ENTER NUMBER OF SEASONAL AR PARAMETERS. 
0 
ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 


u 
ENTER NUMBER OF SEASONAL MA PARAMETERS. 


INITIAL PARAMETER ESTIMATES FOR MARQRT : 
THETA(1) = 0.390425 
THETAS(1) = 0.534397 


DO YOU WANT TO GO AGAIN? 
x 


ENTER LETTER FOR OPTION YOU WANT. 


Ww 
IS YOUR DATA IN FILE FTO2F001? 


sf 
EXECUTION BEGINS... 


ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 
Bsrep NUMBER OF SEASONAL DIFFERENCES. 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 
ENTER NUMBER OF SEASONAL AR PARAMETERS. 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 


H~ oOo OO fF 


ENTER NUMBER OF SEASONAL MA PARAMETERS. 
1 

ENTER LENGTH OF SEASON VIA FORMAT I2. 
12 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 
ENTER NON-SEASONAL MA PARAMETER THETA(1). 
390425 


ENTER SEASONAL MA PARAMETER THETAS(1). 
- 534397 ae 


x 


CONVERGENCE HAS BEEN REACHED IN MAX( 7, 7) ITERATIONS. 


SELECTED OUTPUT FOLLOWS: 


PARAMETER ESTIMATE STANDARD ERROR 
THETA(1) = 0.377152 0.819933D-01 
THETAS(1) = 0.572387 0.780252D-01 


MOVING AVERAGE CONSTANT: THETAO = 0.000291 


a2 


CHI-SQUARE STATISTIC FOR RESIDUAL LACK OF FIT = 29.713571 


DEGREES OF FREEDOM = 38 


PROBABILITY OF EXCEEDING STATISTIC = 0.829438 
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_ 
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DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 
y/ 


EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.009 0.026 -0.129 -0.105 0.078 0.077 -0.036 -0.033 0.103 -0.050 
0.026 -0.020 0.013 90.035 0.067 -0.130 9.052 0.015 -0.093 -0.090 
-0.027 -0.014 0.213 0.009 -0.042 


PARTIAL AUTOCORRELATIONS 


> 


0.009 0.026 -0.130 -0.104 0.088 0.068 -0.072 -0.029 0.151 -0.057 
-0.022 0.018 0.046 0.001 0.053 -0.114 0.069 0.029 -0.129 -0.134 
0.042 -0.014 0.143 -0.025 0.013 


MEAN=0.236005E-02 VARIANCE = 0.139713E-02 


ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / (0,1,1)x(0,1,1)12 model of series g 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

* . 


DO YOU WANT TO GO AGAIN? 
y. Fi 


ENTER LETTER FOR OPTION YOU WANT. 
y = . 
IS YOUR DATA IN FILE FTO2F001? 


EXECUTION BEGINS... 
ENTER LENGTH OF SEASON VIA FORMAT 12. 
ore NUMBER OF NON-SEASONAL DIFFERENCES. 
“nee NUMBER OF SEASONAL DIFFERENCES. 
eaee NUMBER OF NON-SEASONAL AR PARAMETERS. 
“SNTER NUMBER OF SEASONAL AR PARAMETERS. 


ENTER NUMBER OF NON=SEASONAL MA PARAMETERS. 


1 
ENTER NUMBER OF SEASONAL MA PARAMETERS. 
4 


kkk 


INITIAL PARAMETER ESTIMATES FOR MARQRT ; 


PHIE(1) = 0.112699 
THETA(1) = 0.490883 
THETAS(1) = 0.533679 


DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
Ww 
IS YOUR DATA IN FILE FT02F001? 


yy, 
EXECUTION BEGINS... 
=_— NUMBER OF NON-SEASONAL DIFFERENCES. 


ENTER NUMBER OF SEASONAL DIFFERENCES. 
1 


ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 
it 
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ENTER NUMBER OF SEASONAL AR PARAMETERS. 


0 
ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 


1 
ENTER NUMBER OF SEASONAL MA PARAMETERS. 
il 


ENTER LENGTH OF SEASON VIA FORMAT 
12 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, 


2 


ENTER NON-SEASONAL AR PARAMETER PHI(1). 


0.112699 
0.490883 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 


ENTER SEASONAL MA PARAMETER THETAS(1). 


0.533679 


ake 


CONVERGENCE HAS BEEN REACHED IN MAX(10,10) ITERATIONS. 


SELECTED OUTPUT FOLLOWS: 
PARAMETER ESTIMATE 


PHI(1) = 0.146667 
- ° ‘THETA(1) = 0.509816 
THETAS(1) = 0.573388 


MOVING AVERAGE CONSTANT: 


STANDARD ERROR 


0.224464D 00 
0.195934D 00 
0.784896D-01 


THETAO = 0.000248 


AS REQUESTED. 


CHI-SQUARE STATISTIC FOR RESIDUAL LACK OF FIT = 28.302543 


DEGREES OF FREEDOM = 37 


PROBABILITY OF EXCEEDING STATISTIC = 


DO YOU WANT TO PLOT AUTO AND PAUTO 


y 
*EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.013 0.073 -0.100 -0.088 0.077 
0.033 -0.018 0.014 0.027 0.070 
-0.020 -0.023 0.208 0.005 -0.036 


ad 


PARTIAL AUTOCORRELATIONS 
=0.013 0.073 -0.098 -0.096 0.092 
=0.027 0.017 0.047 -0.004 0.053 
0.051 -0.006 0.148 -0.019 0.003 
MEAN=0.253052E-02 VARIANCE 


ENTER TITLE FOR PLOTS. ; 
autos and pautos of residuals / (1 


OF RESIDUALS? 


0.073 -0.031 -0.033 
-0.133 0.051 0.008 


0.081 -0.064 -0.040 
0.118 0.055 0.031 


= 0.139373E-02 


71,1)x(0,1,1)12 model 


YOUR AUTO AND PAUTO PLOTS‘HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
y 
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0.847102 


0.105 -0.054 
-0.092 -0.084 


0.152 -0.050 
-0.123 -0.135 


of series g 
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ENTER LETTER FOR OPTION YOU WANT. 
f 


IS YOUR DATA IN FILE FTO2F001? 
of 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 
y, 


ENTER SEASONAL DATA AS REQUESTED: 


- ENTER 
1 
ENTER 
0 
ENTER 


NUMBER OF SEASONAL DIFFERENCES. 
NUMBER OF SEASONAL AR PARAMETERS. 
NUMBER OF SEASONAL MA PARAMETERS. 


ENTER 
12 


LENGTH OF SEASON VIA FORMAT I2. 


ENTER NONSEASONAL DATA AS REQUESTED: 

' ENTER 
1 
ENTER 


NUMBER OF NON-SEASONAL DIFFERENCES. 


NUMBER OF NON-SEASONAL AR PARAMETERS. 
0 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 
at 
NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, THETAO. 
-000291 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 
=e) Tey 


ENTER SEASONAL MA PARAMETER THETAS(1). 
NO 238 7 


ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT 12. 
30 


ENTER INDEX FOR PLOT ORIGIN VIA FORMAT I3. 
100 


ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT I3. 
ats} 


ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 
0.10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
y 


DO YOU WANT BASIC OUTPUT AT THE TERMINAL? 
Y 


THE LAST 10 WVEC VALUES: 


0.418999E 03 0.461001E 03 0.472000E 03 0.535001E 
0.606000E 03 0.507999E 03 0.461001E 03 0.390001E 
THE 30 FORECAST VALUES: 

0.398539F 03 0.415184E 03 0.395402E 03 0.462697E 
0.469710E 03 0.543456E 03 0.618194E 03 0.626276& 
0.460421E 03 0.405213E 03 0.448016E 03 0.466864E 
0.520593E 03 0.507548E 03 0.528791E 03 0.611991E 
0.705666E 03 0.590927E 03 0.519089E 03 0.456978E 
0.526812E 03 0.502002E 03 0.587782E 03 0.573220E 


£57 


03 
03 


03 
03 
03 
03 
03 
03 


0.622000E 
0.432002E 


0.450971E 
0.524292E 
0.444748E 
0.696357E 
0.505397E 
0.597386E 


03 
03 


03 
03 
ee) 
03 
03 
03 


oem. , 
ew DOV Be 
= 

oi eb Oe ll 


AN 


3eew Of 


ae.0 EG 
elon. Ee 


[a ULAY 


sretlétpul ; 


a 
a 


Po 


THE 30 UPPER 


FORECAST CONFIDENCE LIMITS: 


_ ENTER 


0.399598E 03 0.416254E 03 0.396481E 03 0.463785E 
0.470813E 03 0.544567E 03 0.619311E 03 0.627399E 
0.461557E 03 0.406354E 03 0.449173E 03 0.468030E 
0.521779E 03 0.508743E 03 0.529995E 03 0.613203E 
0.706894E 03 0.592163E 03 0.520332E 03 0.458229E 
0.528090E 03 0.503292E 03 0.589083E 03 0.574532E 
ALPHA FOR THE CONFIDENCE LIMITS IS: 0.100 


® =, 
DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
£ 


Is YOUR DATA IN FILE FTO2F001? 
4 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 
Y 


SEASONAL DATA AS REQUESTED: 


ENTER 
1 

ENTER 
0 

ENTER 
iL 

ENTER 
12 


NUMBER OF SEASONAL DIFFERENCES. 
NUMBER OF SEASONAL AR PARAMETERS. 
NUMBER OF SEASONAL MA PARAMETERS. 


LENGTH OF SEASON VIA FORMAT I2. 


ENTER NONSEASONAL DATA AS REQUESTED: 
ENTER 
1 
ENTER 
1 
ENTER 
1 


NUMBER OF NON-SEASONAL DIFFERENCES. 
NUMBER OF NON=SEASONAL AR PARAMETERS. 


NUMBER OF NON-SEASONAL MA PARAMETERS. 
NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, THETAO. 

-000248 

ENTER NON-SEASONAL AR PARAMETER PHI(1). 
0.146667 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 
0.509816 


ENTER SEASONAL MA PARAMETER THETAS(1). 
0.573388 


ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT I2. 
30 


ENTER INDEX FOR PLOT ORIGIN VIA FORMAT I3. 
100 


ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT I3. 
32 


ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 
0.10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
iy 


DO YOU WANT.BASIC OUTPUT AT THE TERMINAL? 
y 


158 


03 
03 
03 
03 
03 
03 


0.452068E 
0.525422E 
0.445925E 
0.697577E 
0.506664E 
0.598708E 


03 
03 
03 
03 
03 
02 


(CPAgg® 


JAN IHOTY ae 


THE LAST 10 WVEC VALUES: 


0.418999E 
0.606000E 


03 
03 


0.461001E 
0.507999E 


. THE 30 FORECAST VALUES: 


0.398650E, 
0.469348E 
_ 0.460071E 

-0.520096E 
0.704905E 
0.526113E 


03 
03 
03 
03 
03 
03 


THE 30 UPPER 


' 0.399709E 
0.470448E 
0.461199E 
0.521270E 
0.706117E 
0.527371E 


ALPHA FOR 


03 
03 
03 
03 
03 
03 


THE CONFIDENCE LIMITS IS: 0.100 


0.414875E 
0.543045E 
0.404899E 
0.507066E 
0.590336E 
0.501400E 


03 
03 


03 
03 
03 
03 
03 
03 


0.472000E 
0.461001E 


0.395125E 
0.617700E 
0.447739E 
0.528257E 
0.518566E 
0.587074E 


03 
03 


03 
03 
03 
03 
03 
03 


FORECAST CONFIDENCE LIMITS: 


0.415945E 
0.544151E 
0.406032E 
0.508249E 
0.591556E 
0.502668E 


DO YOU WANT TO GO AGAIN? 


¥ 


03 
03 
03 
03 
03 
03 


0.396204E 
0.618812E 
0.448887E 
0.529448E 
0.519792E 
0.588353E 


ENTER LETTER FOR OPTION YOU WANT. 


E 


IS YOUR DATA IN FILE FTO2F001? 


Me 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 


Y 


ENTER SEASONAL DATA AS REQUESTED: 


ENTER 
1 

ENTER 
0 

ENTER 
1 

ENTER 
2 


ENTER 


ENTER 
1 

ENTER 
0 

ENTER 
1 


NUMBER OF SEASONAL DIFFERENCES. 


LENGTH OF SEASON VIA FORMAT I2. 


NONSEASONAL DATA AS REQUESTED: 


NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, 


-000291 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 


~377152 


THETAO. 


ENTER SEASONAL MA PARAMETER THETAS(1). 


572387 


ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT 


40 
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03 
03 
03 
03 
03 
03 


NUMBER OF SEASONAL AR PARAMETERS. 


NUMBER OF SEASONAL MA PARAMETERS. 


NUMBER OF NON-SEASONAL DIFFERENCES. 


0.535001E 
0.390001E 


0.462366E 
0.625753E 
0.466360E 
0.611381E 
0.456512E 
0.572532E 


0.463452E 
0.626870E 
0.467517E 
0.612579E 
0.457745E 
0.573821E 


NUMBER OF NON-SEASONAL AR PARAMETERS. 


NUMBER OF NON=-SEASONAL MA PARAMETERS. 


12. 


03 
03 


03 
03 
03 
03 
03 
03 


03 
03 
03 
03 
03 
03 


0.622000E 
0.432002E 


0.450651E 
0.523897E 
0.444325E 
0.695632E 
0.504960E 
0.596632E 


0.451744E 
0.525020E 
0.445491E 
0.696837E 
0.506207E 
0.597930E 


03 
03 


03 
03 
02 
03 
03 
03 


03 
03 
03 
03. 
03 
03 
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o> TRADE 


ENTER INDEX FOR PLOT ORIGIN VIA FORMAT 13. 

075 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT 13. 

120 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 
0.10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
Y 


DO YOU WANT BASIC OUTPUT AT THE. TERMINAL? 


YOU WANT TO GO AGAIN? 


eae ke 


ENTER LETTER FOR OPTION YOU WANT. 
d 
IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


,; ENTER ORDER OF SEASONAL DIFFERENCING. 

1 

ENTER LENGTH OF SEASONAL PERIOD VIA I[2. 
12 


ENTER NUMBER OF NONSEASONAL DIFFERENCES. 
1 


DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


n 
DO YOU WANT TO GO AGAIN? 
y; 


ENTER LETTER FOR OPTION YOU WANT. 
e 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

alter file £t02f£001 pl Inserg £t02f001 pl 
alter file £t03f£001 pl file £t02f001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

Mf 


ENTER LETTER FOR OPTION YOU WANT. 
x 


IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 


y/ 
ENTER LENGTH OF SEASON VIA FORMAT 12. 
12 


ENTER NUMBER OF SEASONAL AR PARAMETERS. 
0 


ENTER NUMBER OF SEASONAL MA PARAMETERS. 
1 


ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 
0 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 
1 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 
-377152 


ENTER SEASONAL MA PARAMETER THETAS(1). 
-572387 
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FOR THE GIVEN INPUT PARAMETERS AND MODEL 


THE RESIDUAL SUM OF SQUARES IS: 0.18192935D 00 


DO YOU WANT TO TEST DIFFERENT PARAMETER VALUES? 
n ; 


DO YOU WANT TO TEST A DIFFERENT MODEL? 
y 


IS YOUR SERIES SEASONAL? 


VENTER LENGTH OF SEASON VIA FORMAT I2. 
Seerek NUMBER OF SEASONAL AR PARAMETERS. 
Seer en NONSER OF SEASONAL MA PARAMETERS. 
Parer NUMBER OF NON-SEASONAL AR PARAMETERS. 


1 
_ ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 
Pi 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 


ENTER NON-SEASONAL AR PARAMETER PHI(1). 

- 146667 

ENTER NON-SEASONAL MA PARAMETER: THETA(1). 
-509816 

ENTER SEASONAL MA PARAMETER THETAS(1). 
573388 


* 


FOR THE GIVEN INPUT PARAMETERS AND MODEL 


THE RESIDUAL SUM OF SQUARES IS: 0.18164378D 00 


DO YOU WANT TO TEST DIFFERENT PARAMETER VALUES? 
n 


DO YOU WANT TO TEST A DIFFERENT MODEL? 
n 


DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 
R; 
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PARTIAL AUTOCORRELATIONS 


163 


14292.0 
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PROGRAM LISTINGS 


This Appendix contains program listings of all programs 
resident in the Time Series Editor. It does not include 
listings of subroutines used from other libraries, such 


as IMSL routines. 
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APPENDIX E 
LISTING OF DATA SETS RESIDING 
IN THE TIME SERIES EDITOR 
This appendix contains listings of the data sets now 

residing in the Time Series Editor. The data consists of 
Box and Peaking? [Ref. 4] time series C [Ref. 4, p. 528] 
and time series G [Ref. 4, p. 531], as well as a data set 
containing monthly Monterey, California, rainfall data 


from January 1931 through December 1976 . 
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